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PREPACE 


This  report  describes  the  work  performed  at  Physical  Sciences  Inc.  (PSI)  as  a 
Phase  I  Small  Business  Innovative  Research  (SBIR)  program  sponsored  by  the 
Office  of  Naval  Research:  The  general  goal  of  the  program  was  to  demonstrate 
that  the  fluid  flow  measurement  technique  known  as  the  Vorticity  Optical 
Probe  (VOP)  could  be  used  to  measure  the  vorticity  vectors  within  a  volume  of 
a  flow  field  at  many  distinguishable  locations  simultaneously.  The  spatial 
resolution  of  the  measurements  would  be  small  enough  for  them  to  be  sensitive 
to  the  fine-scale  fluctuations  of  the  flow,  yet  the  overall  sampling  volume 
large  enough  to  enable  study  of  the  interactions  between  fine-scale  features 
and  large-scale  structures.  The  specific  objectives  of  the  program  were  to 
demonstrate  that:  1)  The  VOP  provides  sufficient  information  to  allow  rapid 
measurement  of  all  three  components  of  the  vorticity  vector  at  a  single  well- 
defined  point  within  the  flow  field;  2)  The  vorticity  vectors  at  several 
spatially  distinct  but  unknown  positions  can  be  measured  concurrently;  and  3) 
A  method  could  be  developed  for  determining  where  in  the  sampling  volume  the 
vorticity  is  being  measured. 


These  three  specific  goals  have  been  met  using  simultations  of  both  the  VOP 
and  rotational  flow.  A  data  analysis  technique  has  been  developed  which  demon¬ 
strates  the  three-component  measurement  capability  and  indicates  the  inherent 
limits  on  its  precision.  The  capability  to  obtain  vorticity  data  from  many 
points  simultaneously  has  been  demonstrated,  and  it  was  found  that  there  is  a 
correlation  between  the  accuracy  of  the  vorticity  measurement  and  the  preci¬ 
sion  with  which  the  location  of  the  measurement  within  the  large  sampled  vol¬ 
ume  is  known.  In  this  Phase  I  program  we  have  shown,  as  planned,  that  one  of 
the  three  position  coordinates  needed  to  pinpoint  the  measurement  location  can 
be  easily  determined,  and  that  the  other  two  can  in  principle  be  found  by  a 
conceptually  straightforward  extension  of  the  one-coordinate  technique.  We 
have  also  conceived  some  other  position  measurement  methods  that  we  were 


V  .. 

V-,' 


WE 


ss-a 


m 


Ml 


■\v.yi 


'.'.'I 


>>./ 


unable  to  pursue  in  Phase  I  but  believe  will  ultimately  be  relatively  easy  to 
incorporate  into  a  volumetric  VOP  system. 

The  remainder  of  this  report  describes  the  details  of  these  accomplishments. 
Section  1  is  primarily  an  introductory  section  which  includes  the  background 
information  needed  for  comprehension  of  subsequent  material.  Section  2  pre¬ 
sents  the  Phase  I  work  and  the  specific  tasks  which  were  performed  to  meet  the 
objectives  listed  above.  The  report  concludes  with  a  brief  discussion  of  the 
hurdles  which  we  would  intend  to  overcome  in  a  Phase  II  program  leading  to  com¬ 
mercialization  of  the  VOP. 
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1.  VOP  FUNDAMENTALS 


1.1  INTRODUCTION 

Experimental  fluid  mechanics,  can  be  roughly  divided  into  two  areas;  (1)  Flow 
visualization  (volumetric)  experiments,  wherein  large  volumes  of  a  flow  field 
are  marked  by  dyes,  smoke  are  other  imhomogenei ties ,  photographed,  and  inter¬ 
preted  by  eye,  and  (2)  Probes  of  properties  at  discrete  points  (principally 
anemometry).  Flow  visualization  has  been  successful  in  educing  gross  features 
of  coherent  structures,  but  yields  little  quantitative  information.  In  con¬ 
trast,  anemometry  provides  quantitative  information  at  many  scales  and,  with 
clever  averaging  procedures,  can  describe  some  aspects  of  characteristic  large 
scale  flow  structures  and  their  effects  on  the  surrounding  flow  field. 

Several  efforts  have  been  made  to  derive  quantitative  information  from  sophis¬ 
ticated  flow  visualization  apparatus,  but  the  data  analysis  requirements 
overwhelm  even  the  largest  computers.*  Unfortunately,  even  the  most  exacting 
measurements  have  been  unable  to  achieve  sufficient  local  detail  over  large 
enough  flow  volumes  to  show  the  interaction,  if  any,  of  large  and  small  scale 
structure.  Conditional  averaging  procedures?  for  continuous  measurements  at 
discrete  points  necessarily  filter  out  information  in  all  but  a  narrow  band  of 
large  length  scales. 

Quantitative  volumetric  measurements  of  flow  properties  over  a  broad  range  of 
length  scales  are  sorely  needed  to  advance  our  understanding  of  complex  flows 
and  ultimately  develop  the  ability  to  alter  and  control  them.  Specifically, 
techniques  are  required  that  offer  sufficient  spatial  resolution  to  distin¬ 
guish  fine  scale  motions  over  a  sufficiently  large  volume  to  determine  how  the 
fine-scale  features  interact  to  form  large  scale  structures,  or  how  the  fine 
scales  are  formed  by  interactions  of  the  large  scale  structures  within  the 
mean  flow. 


The  essential  distinction  between  the  two  areas  of  experimental  fluid  dynam¬ 
ics,  flow  visualization  and  discrete  fixed  probes,  is  that  the  first  is  an 
attempt  to  follow  the  motion  of  many  (marked)  fluid  elements,  and  the  second 
measures  some  property  of  the  flow  as  it  convects  past  fixed  points; 

Lagrangian  and  Eulerian  respectively.  Few  experimental  techniques  fall  into 
the  middle  ground  of  this  global  division.  As  described  below,  the  volumetric 
Vorticity  Optical  Probe  would  be  such  a  technique.  The  individual  probe  parti 
cles  are  passive  co-moving  markers  of  the  flow  field,  but  at  the  same  time, 
they  are  active  probes  of  the  local  vorticity  vector  and  offer  a  fine  intrin¬ 
sic  spatial  resolution  of  -25  pm.  The  complementary  nature  of  the  data  thus 
obtained,  with  appropriate  analysis,  can  potentially  provide  more  information 
about  the  structure  of  the  sampled  flow  than  comparable  volumes  of  vorticity 
or  velocity  data  alone. 

1.2  THE  VOP 

Appendix  A  is  a  copy  of  a  published  paper  which  provides  the  essential  details 
of  the  technique  recently  developed  by  Frish  and  Webb^  for  direct  optical 
measurements  of  vorticity.  The  technique,  known  as  the  Vorticity  Optical 
Probe  or  VOP,  is  based  on  the  propensity  of  spherical  particles  suspended  in  a 
flowing  fluid  to  rotate  with  an  angular  velocity  which  is  equal  to  one-half 
the  local  vorticity.  Since  the  20  -  50  ym  diameter  spheres  developed  for  this 
purpose  are  smaller  than  the  smallest  scale  in  liquid  turbulent  flows  (the 
Kolmogorov  microscale)  and  have  small  masses  and  moments  of  inertia,  they  have 
sufficiently  rapid  response  times  to  react  to  the  fastest  fluctuations  of 
vorticity  and  velocity  in  the  flow. 

To  measure  the  angular  velocity  of  such  small  particles,  specially  manufac¬ 
tured  polymethylmethacrylate  (PMMA  or  Plexiglas)  spheres  which  contain  one  or 
more  embedded  lead  carbonate  crystals  are  used.  The  crystals  are  planar, 
highly  reflective  hexagons,  typically  16  ym  wide  and  0.06  ym  thick.  When 
suspended  in  a  liquid  which  matches  the  index  of  refraction  of  the  PMMA,  the 
spheres  are  invisible.  A  beam  of  light,  such  as  from  a  laser,  impinging  on 


such  a  sphere  is  reflected  by  the  lead  carbonate  crystal  in  a  direction  deter¬ 
mined  by  the  relative  orientation  of  the  crystal  surface  with  respect  to  the 
direction  of  the  light  beam's  propagation.  When  the  sphere  rotates  due  to  the 
presence  of  vorticity,  the  direction  of  the  reflected  beam  changes  with  time. 
By  placing  an  optical  detector  within  the  path  of  the  reflected  beam  and 
measuring  temporal  changes  in  the  direction  towards  which  it  is  reflected,  as 
indicated  in  Figure  1,  the  angular  velocity  of  the  sphere  and  thus  the  vorti- 
city  of  the  fluid  is  deduced. 

1.2.1  VOP  Geometry 

The  coordinate  system  shown  in  Figure  2  illustrates  the  operation  of  the  VOP. 

A  light  beam,  incident  along  the  x-axis,  approaches  the  origin  where  a 
reflective  particle  is  located.  The  mirror  orientation,  defined  as  the 
direction  of  its  normal  vector,  n,  is  given  by  the  two  angles  9  and  <|>.  The 
reflections  emerge  in  the  n-x  plane  at  an  angle  defined  in  these  coordinates 
as  (9',  V ) .  The  vorticity  of  the  particle  (w*,  Wy,  u^)  is  determined  by 
measuring  and  •  To  relate  these  measured  values  to  w,  we  first  find  the 

relationship  between  (9,  <t>)  and  (9', 


Figure  1. 


Illustration  of  trajectory  along  a  spherical  surface  of 
reflection  from  a  mirror  rotating  with  vorticity  w. 


INCIDENT  LIGHT 
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Figure  2.  Coordinate  system  for  analysis  of  reflection  trajectory. 
It  is  easiest  to  work  entirely  in  cartesian  coordinates.  Thus 

x  =  sin0cos<f> 
y  =  sin9sin<|> 
z  =  cos©. 

The  x-axis  and  n  define  a  plane  which  forms  an  angle  C  with  the  x-z  plane 
which  contains  both  the  incident  and  reflected  beams.  If  y  is  the  angle 
between  the  x-axis  and  n  on  this  plane,  then  2y  is  the  direction  of  the 
reflection.  Trigonometry  yields 
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y  =  sinysinC 
z  =  sinycosC 


Equating  Eqs.  (la-c)  with  (2a-c)  yields 

cosy  =  sin©cos4> 

-sinysinC  =  sin0sin<J> 
sinycosC  =  cos0, 

and  solving  for  0  and  <|>  in  terms  of  y  and  C  gives 

0  =  cos_^(sinycosC) 

<f>  =  cos-1(-cosy/  Jl-sin^ycos^C) 

Thence 

0'  =  cos"l(sin2ycosC) 

<)>'  =  cos“^(cos2y/  Jl-sin^ycos^C) 

Solving  now  for  y  and  C  in  terms  of  0  and  and  substituting  back  into  Eqs 
(5a  and  5b)  gives 

0'  =  cos-l(sin20cos<j>) 

<f>'  =  cos-l((sin2ecos2.-cos20)/JI—i^;] 
which  are  the  desired  relationships. 
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Now,  if  9  and  $  vary  with  time,  the  rates  of  change  of  9'  and  <J>'  are  found  to 
be 


and 


d9'/dt  =  [1/  Jl-cos2$sin220] 

x  [  (sin20sin<^)d4>/dt-(2cos20cos<j>)d9/dt)  J 


d^'/dt  =  [l/(l-cos24>sin220) 

x  [  (2sin29-sin220cos24>)d4>/dt-(sin29sin24>)d0/dt  ] 


The  rate  of  change  of  n  is  related  to  the  vorticity  components  by^ 

d9/dt  =  WyC0s4>  -  o^sin^ 

d4>/dt  =  -  co t9( (0x003  41  +  <0ysin4>) 


(7a) 


(7b) 


(8a) 

(8b) 


When  an  optical  detector  records  the  reflection  passing  through  the  orienta¬ 
tion  (9',  4>' )  relative  to  the  incident  beam,  a  unique  relationship  between  the 
components  of  the  vorticity  vector  and  the  local  values  of  d9/dt  and  d4>/dt 
(i.e.,  the  local  reflection  trajectory)  is  established.  Vorticity  vector 
components  can  therefore  be  calculated  by  measuring  d9'/dt  and  d4>'/dt  at  that 
orientation.  For  example,  at  9  =  n/2,  and  $  =  n/4; 


and 


so  that 


and 


d0/dt  =  (wy  -  wz)/2'f2 
d4»/dt  =  0^/2  , 

d9'/dt  =  (toy  -  tox)/,T?  , 
d$'/dt  =  2wz  . 


(9a) 


(9b) 


(10a) 


(10b) 
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1.2.2  Trajectory  Simulations 


In  the  general  case  of  a  large  solid  angle  detector  surface,  Eqs.  (7)  and  (8) 
do  not  reduce  to  the  simple  form  of  Eqs.  (9)  and  (10)  since  0(n)  and  <J>(n)  can 
take  on  large  ranges  of  values  over  which  the  reflections  will  strike  the 
detector.  The  passage  of  the  reflected  beam  across  the  detector  surface  is 
referred  to  as  a  reflection  trajectory,  or  just  trajectory.  For  any  given 
vorticity  vector,  there  is  a  corresponding  family  of  trajectories,  each  member 
of  which  has  a  different  initial  mirror  normal  vector,  n,  and  thus  traverses  a 
different  path  on  the  detector.  Figure  3(a-d)  shows  several  simulated  trajec¬ 
tory  families  for  selected  vorticity  vectors.  The  trajectories  traverse 
across  a  flat  1  sr  solid  angle  surface.  The  angular  coordinates  of  the  surface 

(detector)  center  are  0  =  <f>  =  £  or  0'  =  Figure  3(a)  shows 

possible  trajectories  for  pure  <0%,  i.e,  =  Wy  =  0.  Figure  3(b)  adds  a  com¬ 
ponent  perpendicular  to  and  parallel  to  the  "principal  reflection  plane", 
defined  as  that  plane  which  would  reflect  the  incident  beam  to  the  detector 
center.  The  effect  of  this  additional  vorticity  component  is  mainly  to 
incline  the  trajectories  across  the  surface.  Figure  3(c)  includes  yet  another 
component,  now  perpendicular  to  this  plane  and  thus  also  perpendicular  to  the 
two  components  of  Figure  3(b).  Note  that  the  differences  between  Figures  3(b) 
and  3(c)  appear  almost  entirely  as  subtle  changes  of  trajectory  curvatures. 

If  this  third  component  should  become  dominant  as  in  Figure  3(d),  reflections 
sweep  about  a  coordinate  singularity  passing  through  the  detector.  Clearly, 
to  analyze  trajectories  to  yield  all  three  components  of  vorticity  using  a 
single  detector  surface,  it  is  essential  to  not  simply  record  the  rates  at 
which  reflections  pass  through  the  detector's  center  point,  but  to  also  deter¬ 
mine  precise  the  local  curvature  of  the  trajectories. 


1 . 3  PREVIOUS  WORK 
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Prior  to  the  current  program,  all  work,  on  the  VOP  was  performed  by  the  two 
authors  of  this  report  during  their  respective  graduate  careers  at  Cornell 
University.  Originally,  the  VOP  was  tested  using  the  optical  system  illus¬ 
trated  in  Figure  4.  To  measure  a  single  component  of  the  vorticity  vector 
over  a  small  volume,  which  can  be  considered  a  single  point,  a  laser  beam 
enters  the  flow  channel  at  a  45  degree  angle  through  a  prism  used  to  prevent 
refractior.  of  the  beam.  The  flow  channel  is  filled  with  a  small  concentration 
of  randomly  oriented  VOP  particles  suspended  in  a  liquid  which  matches  the 
refractive  index  of  the  PMMA  spheres.  Usable  vorticity  measurements  are 
generated  by  those  particles  which  have  mirrors  oriented  so  as  to  reflect  the 
laser  beam  at  a  nearly  perpendicular  angle  within  the  plane  defined  by  the 
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Figure  4.  Original  VOP  optical  detection  system. 
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incident  beam  and  the  detecting  optics.  These  particles  are  imaged  onto  slit 
3  which,  together  with  the  incident  beam,  define  the  sampled  volume.  The 
rotational  velocity  of  the  reflected  light  from  each  properly  oriented  mirror 
is  measured  by  observing  its  transit  time  across  a  small  angle  defined  by  a 
pair  of  closely  spaced  slits,  shown  as  slits  1  and  2.  Light  passing  through 
these  slits  is  converted  to  electrical  impulses  by  photomultiplier  tubes 
(PMTs).  The  time  between  reception  of  sequential  pulses  from  the  two  PMTs  is 
a  measure  of  the  vorticity. 

This  system  was  tested  in  a  laminar  Poiseulle  flow  (see  Appendix  A)  and  used 
to  measure  spanwise  (a^)  vorticity  fluctuations  and  profiles  in  a  transitional 
boundary  layer. ^  This  system  was  restricted,  however,  to  measurement  of  only 
the  (dominant)  spanwise  component  of  vorticity  (u^)  in  flows  having  large 
vorticity/velocity  ratios.  Also,  the  working  fluid's  refractive  index  was 
required  to  match  PMMA  (1.49).'  Organic  liquids,  specifically  dibutylphthalate 
and  p-cymene,  were  used  for  that  purpose. 

The  second  generation  optical  configuration, ^  shown  in  Figure  5,  addressed  the 
restriction  on  the  vorticity/velocity  ratio,  but  still  required  exotic  fluids 
for  index  matching.  The  slits  and  photo-multipliers  of  the  original  system 
were  replaced  by  a  two-dimensional  position-sensing  photodiode  enabling  the 
simultaneous  measurement  of  spanwise  (o^)  and  streamwise  (ws)  vorticity  compon¬ 
ents.  The  sampled  volume  was  imaged  onto  the  pinhole  shown  in  the  figure,  and 
reflections  traversed  the  photodiode.  The  photodiode  produced  analog  signals 
corresponding  to  the  (x,y)  position  of  the  reflection  on  the  detector  surface, 
providing  not  only  a  rate  of  passage,  but  the  inclination  of  the  trajectory. 
Near  the  detector  center  this  information  gives  accurate  estimates  of  two  com¬ 
ponents  of  vorticity  without  requiring  a  minimum  angular  displacement,  thus 
eliminating  vorticity/velocity  ratio  restrictions.  With  an  improved  optical 
system,  the  second  generation  VOP  could  potentially  yield  all  three  vorticity 
components . 


This  second  generation  configuration  was  used  to  obtain  simultaneous  single¬ 
point  measurements  of  and  «s  in  a  turbulent  boundary  layer.  Considerable 
streamwise  vorticity  was  observed.  The  limitations  of  this  system  were  lack 
of  water  compatibility,  and  difficulty  in  discriminating  between  single  (good) 
and  multiple  (bad)  reflection  events  thus  requiring  elaborate  off-line  data 
analysis.  Multi-point  spatial  vorticity  correlations  made  with  two  such 
systems  imaging  adjacent  points  were  confounded  by  tempqrally  uncorrelated 
sampling,  resulting  in  poor  statistics. 


DUAL  AXIS 
POSITION-SENSING 
PHOTODIODE 
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Figure  5.  Second  generation  optical  system  for  two-component 
vorticity  measurement  (from  ref.  6). 
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2.  PHASE  I  PROGRAM 


As  we  had  proposed,  the  Phase  I  program  was  divided  into  three  distinct  tasks. 
Task  1  demonstrated  that  all  three  components  of  the  vorticity  vector  can  be 
deduced  from  a  single  reflection  trajectory  across  a  single  detector  surface. 
Task  2  showed  that  many  trajectories  can  be  captured  simultaneously,  and  be 
interpreted  for  data  reduction.  Task  3  was  a  demonstration  of  the  feasibility 
of  establishing  the  origins  from  which  individual  reflection  emanate  within 
the  sampled  volume.  All  of  these  tasks  have  been  completed  successfully, 
thereby  demonstrating  the  possiblity  of  assembling  a  system  capable  of 
volumetric  vorticity  measurement. 

2.1  TASK  1  -  THREE  COMPONENT  VORTICITY  MEASUREMENT 

As  discussed  above,  previous  work  on  the  VOP  was  restricted  to  measurements  of 
one  or  two  components  of  the  vorticity  vector,  at  one  or  two  fixed  positions 
in  the  flow.  The  simulated  reflection  trajectories  of  Figure  3  indicate, 
however,  that  careful  analysis  of  the  actual  trajectories  can  yield  all  three 
vorticity  components.  In  contrast  to  the  concept  of  using  two  or  more  widely 
separated  detectors  as  we  had  proposed  for  this  measurement,  the  simulations 
showed  that  unambiguous  vorticity  values  could  be  determined  by  making  many 
measurements  of  dQ' /dt  and  d$'/dt  at  modest  angular  separations  along  the 
surface  of  a  single  large-solid-angle  detector.  Therefore,  for  this  task,  a 
single  photographic  plate  was  used  as  the  detector. 

2.1.1  Apparatus 

Verification  of  three  component  vorticity  measurement,  its  precision,  and  its 
sensitivity  to  various  component  mixtures  and  detection  geometries,  was  most 
easily  accomplished  with  a  VOP  simulation  using  a  plane  mirror,  as  in  the  appa¬ 
ratus  shown  in  Figure  6a.  The  motor  acted  as  an  adjustable  vorticity  vector 
(magnitude  and  direction)  and  a  plane  mirror  affixed  to  the  shaft  end  served 
as  the  mock  VOP  mirror.  A  chopped  Argon  laser  beam  (514.4  nm)  was  incident 
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along  the  positive  x-axis.  Photographic  film  was  was  used  as  detector.  The 
film  plane  was  parallel  to  the  x-z  plane,  at  any  y  value  desired.  The  detec¬ 
tion  geometry  was  fixed  in  this  way  throughout  Task  1. 


B 


The  chopped  incident  beam  provided  a  time  base  so  that,  as  the  motor  shaft 
turned,  the  reflections  produced  a  sequence  of  spots  or  streaks  on  the  film  at 
known  time  intervals.  The  two-angle  variable  motor  mount  with  variable  speed 
motor  allowed  any  mixture  of  the  angular  velocity  components,  and  thus  simula¬ 
ted  vorticity  components  w*,  Wy,  relative  to  the  detection  geometry  while 
maintaining  the  mock  VOP  mirror  precisely  at  the  origin.  Because  the  rotation 
axis  can  have  any  orientation  relative  to  the  incident  beam,  the  orientation 
of  mirror  normal  was  made  adjustable  in  order  to  insure  that  the  reflections 
cross  the  film  plane.  Figure  6b  shows  details  of  the  variable  pitch  mirror 
mount.  The  various  orientations  produce  the  individual  members  of  trajectory 
families  as  described  in  Sub-section  1.2.2.  Note  that  the  measured  vorticity 
must  be  independent  of  the  orientation  of  the  mirror  normal. 
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a)  Task  1  apparatus  for  VOP  simulation. 
Figure  6.  The  rotating  mirror. 
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b)  Adjustable  pitch  mirror  mount. 

Figure  6.  Continued 

The  choice  of  photographic  film  as  a  detector  was  the  logical  one,  since  it 
provided  a  plane  surface  of  large  solid  angle  without  imaging,  and  a  conven¬ 
ient  digitizable  record  for  data  analysis.  The  photographic  film  was 
4x5  in.  Kodax  TMax  400  black  and  white,  mounted  in  a  double  sided  holder 
such  that  the  film  center  corresponded  to  (0,  Yf^m,  0)  in  the  geometry  of 
Figure  6a.  In  all  Task  1  experiments,  YfQm  was  selected  to  be  -10  cm  so  that 
the  solid  angle  subtended  by  the  film  plane  was  -1  sr.  The  film  was  developed 
in  our  own  darkroom,  and  reflection  positions  were  electronically  digitized 
using  a  digitizer  pad  coupled  to  an  IBM-PC. 

The  entire  apparatus  of  Figure  6a  was  placed  in  a  light-tight  enclosure.  The 
film  holder  was  mounted  so  that  it  could  be  remotely  opened  and  closed.  The 
laser  entered  the  enclosure  through  a  variable  speed  shutter.  For  each  trial, 
the  shutter  speed  was  set  to  one  motor  rotation  period.  Thus  the  reflection 
crossed  the  film  plane  once  and  only  once.  With  the  essential  V0P  geometry 
abstracted  in  this  way,  this  configuration  provided  a  controlled  test  of  our 
data  analysis  procedure  and  established  the  inherent  precision  of  the  method. 


2.1.2  Data  Analysis 


Recorded  reflection  trajectories  consisted  of  a  sequence  of  N  spots  on  the 
film.  The  film  negatives  were  taped  onto  the  digitizer  pad,  and  the  spot 
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centroids  or  edges  were  used  for  digitization.  The  position  of  the  if^  spot 
in  the  sequence  (X^,  Yf^m,  Z^)  was  converted  to  reflection  angular  coordi¬ 
nates  (6'i,  $'i)  during  data  reduction,  j  and  j  were  computed  from  spot 


pairs: 


fd©')  0i+l  9i-l  (dr )  *i  +  l  *1-1  u  Ak  .  tU 

Ji  =  — m —  ’  lit  Ji  =  — m — ’  vhere  15  the  chopped 


interval.  Digitization  errors  produce  uncertainties  in  the  computed  deriva¬ 
tives  of  ~1  s_l. 

At  least  three  equations  are  required  to  deduce  a  unique  vorticity  vector. 

('d©' 1  (d*'! 

The  expressions  (7a)  and  (7b)  for  ^--J^and  can  be  written 


(gf-).=  aiWx  +  bi“*  +  ciwz 

(a|-).=  di“x  +  ei“y  +  fi“z 


(11a) 


(lib) 


where  the  coefficients  a-f  are  functions  of  ©'  and  r  which  are  determined 
entirely  by  the  reflection  geometry.  Each  trajectory  produces  2(N-2)  such 
equations.  To  take  full  advantage  of  this  information  and  produce  the  best 
“x.yjZ  values  which  describe  a  trajectory,  a  linear  least  squares  multiple 
parameter  fitting  algorithm  was  written.  The  data  reduction  program  is 
included  as  Appendix  B. 
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2.1.3  Measurements  and  Results 


A  series  of  photographs  were  made  with  several  «  orientations  and  magnitudes. 
The  test  cases  were: 

1.  ft)*,  toy,  *  0 

2.  w*,  toy  *  0;  wz  =  0 

3.  Wj,  *  0;  «*,  «y  =  0 

4.  0%,  to x  =  0 ;  Wy  t  0 

5.  Wy  *  0;  to*  =  0 


g 
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These  test  cases  demonstrate  the  capability  to  measure  any  vorticity  vector 
orientation.  They  included  the  difficult  example  wherein  w*  =  -«y  and  reflec¬ 
tions  sweep  around  a  coordinate  singularity  passing  through  the  film  plane,  as 
in  Figure  3d.  Of  the  data  collected,  two  photographs,  Figures  7  and  8,  have 
been  selected  as  representative  of  best  and  worst  case  precision. 

Figure  7  is  actually  two  trajectories  on  the  same  piece  of  film  made  sequentia 
lly  with  differing  mirror  pitches.  This  photograph  is  case  1;  a^,  <0y,  t  0. 
For  both  pitches,  the  motor  speed  was  10  Hz,  and  the  motor  shaft  orientation 
0S  =  20°,  <frs  =  140°.  For  each  pitch,  the  data  reduction  program  generates  the 
vorticity  component  values  and  their  uncertainties,  via  the  least  square  fit 
to  the  trajectory.  These  u>  values  are  then  converted  back  into  speed,  0S,  and 
♦s  for  comparison.  Table  1  compares  the  vorticity  components  calculated  by 
this  algorithm  with  the  values  expected  from  knowledge  of  the  motor  speed  and 
orientation.  Clearly  the  agreement  is  excellent. 


PITCH  1 


PITCH  2 


Figure  7.  Photograph  of  two  trajectories  at  different  mirror  pitches. 
©-  =  20°,  =  144°. 


TABLE  1.* 


Vorticity  component  calculated  from  trajectories  in  Figure  7. 


«x(s_1) 

Wy(S-1) 

Speed,  0S,  *s 

Pitch  1 

-110.5  ±2 

34.7  ±4 

-29.1  ±4 

9.05  Hz,  22.3°,  140° 

Pitch  2 

-123.6  ±2 

34.6  ±4 

-26.8  +4 

10.45  Hz,  19.5°,  142.2 

Expected 

Values 

-118.0  ±6 

34.7  ±4 

-25.3  ±4 

10.00  Hz,  20.0°,  144° 

*Chopping  rate  600  Hz 

Figure  8  is  case  (2),  (^  =  0,  S  -«y.  Extreme  trajectory  curvature  is  appa¬ 
rent.  Again,  two  different  mirror  pitches  are  included  in  a  single  photo,  but 
different  motor  speeds  were  used.  The  results  are  given  in  Table  2.  Note 
the  vorticity  measurement  in  this  case  establishes  the  orientation  of  the 
motor  shaft  relative  to  the  detection  geometry  (9S,  <|>s)  with  greater  precision 
than  direct  measurement,  due  to  the  large  number  of  points  and  large  trajec¬ 
tory  angular  size. 

2.1.4  Error  Analysis 

The  results  of  Task  1  clearly  demonstrate  three  component  vorticity  measure¬ 
ment  capability  under  ideal  conditions.  The  principal  source  of  error  in 
these  measurements  is  digitization  error.  The  precision  with  which  the 
position  of  the  centroid  of  each  spot  or  streak  on  the  film  is  established 
determines  the  fundamental  precision  of  the  vorticity  measurement,  which  then 
scales  with  both  the  number  of  points  and  trajectory  angular  size.  The  curva¬ 
ture  of  the  trajectory,  which  must  be  measured  to  determine  the  ( co^  -  wyj/'IT 
component,  is  most  highly  sensitive  to  this  source  of  error. 


p; 

('  r) 

t-itch  2 
(Inner) 

Expected 

Values 


Ms-1) 

tH 

1 

to 

sT 

Wy(S_1) 

2.4  +0.3 

73.8  ±0.7 

-77.7  ±0.4 

3.1  ±0.3 

129.6  ±0.6 

-127.2  ±0.6 

1.  0.0  ±2 

88.9  ±6 

-  88.9  ±6 

2.  0.0  ±2 

133.3  ±8 

-133.3  ±8 

8.8  Hz,  88.8°,  135.4' 


10  Hz 


15  Hz 


90.0°,  135° 


*Chopping  rate  600  Hz 

**Motor  speed  drifts  from  nominal  values  during  trails.  0's,  4>'s  independent 
of  motor  speed. 


The  uncertainties  of  the  vorticity  components  are  typically: 


Ato2  <  2 X  of  |<o|,  and 

Af*^,  A«y  ~  2AW2 


The  X^v  goodness  of  fit  statistics  were  generally  <  1.  This  level  of  preci¬ 
sion  or  better  was  observed  in  all  five  tests  cases.  Precision  suffers  with 
diminishing  trajectory  angular  size,  but  could  be  improved  with  more  accurate 
digitization  methods.  Misalignment  of  the  incident  beam  or  film  plan 
introduced  small  additional  systematic  errors.  Motor  speed  fluctuations  also 
enter  from  photo  to  photo,  but  do  not  affect  precise  determination  of  0S  and 
<PS,  which  are  independent  of  motor  speed  and  serve  as  a  basis  of  comparison. 

The  dynamic  range  in  |w|  is  -10.  This  is  determined  by  the  ratio  of  the 
maximum  acceptable  angular  separation  of  "spots"  to  the  minimum  resolvable 
angle  between  them.  This  limitation  is  imposed  by  the  use  of  photographic 
film  and  resolving  spots  by  eye. 

Having  developed  the  analysis  procedure  and  demonstrated  the  inherent  preci¬ 
sion  of  three-component  vorticity  measurement,  we  proceeded  to  study  multiple 
simultaneous  trajectories  from  actual  VOP  mirrors  in  Task  2. 

2.2  TASK  2  -  MULTIPLE  TRAJECTORIES 

The  remainder  of  the  Phase  I  program  was  concerned  with  the  development  of  vol¬ 
umetric  vorticity  measurement  capability.  In  Task  2,  we  simultaneously  recor¬ 
ded  the  trajectories  of  reflections  from  many  lead  carbonate  (VOP)  mirrors  in 
a  large  volume.  VOP  mirrors  embedded  in  a  solid  PMMA  rod  undergoing  steady 
solid  body  rotation  all  have  the  same  angular  velocity  vector.  However,  since 
the  mirrors  are  randomly  oriented,  no  two  reflections  will  necessarily  follow 
the  same  trajectory.  Hence,  this  task  provided  a  stringent  test  of  the  capa¬ 
bility  for  three-component  vorticity  measurement  at  many  points  within  the 
sampled  volume. 


A  mirror  bearing  PMMA  rod  was  manufactured  by  mixing  a  small  quantity  of 
Nacromer  ZTX-B  pearl  pigment  (containing  lead  carbonate  mirrors),  supplied  by 
the  Mearl  Corporation,  with  the  methylmethacrylate  monomer.  Bulk  polymeri¬ 
zation  was  carried  out  in  test  tubes  using  standard  procedures. ^  Several  rods 
with  different  mirror  concentrations  were  produced.  A  concentration  of  ~1(P 
mirrors  per  cm^  was  found  to  be  most  suitable  for  Task  2.  The  rod  was  machined 
to  7  mm  diameter  and  mounted  on  the  motor  in  place  of  the  single  mirror  of 
Task  1. 

As  shown  in  Figure  9,  to  minimize  both  scattering  at  the  rod  surface  and 
refractive  deviation  of  reflections,  the  mirror-bearing  rod  was  suspended  in  a 
refractive  index-matching  fluid  chamber  that  was  attached  to  a  large  Plexiglas 
block. 
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Figure  9.  Task  2  apparatus  for  multiple  trajectory  photographs. 
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The  film  holder  was  mounted  as  close  as  possible  to  the  block  face.  The  index 
matching  fluid  was  dibutylphthalate  (DBT).  The  front  and  rear  Plexiglas  cham¬ 
ber  windows  allow  the  passage  of  the  incident  beam.  The  laser  power  was  ~0.5 
vatts,  and  the  beam  diameter  was  4  mm  giving  an  illuminated  volume  of  0.1  cm^. 
In  all  other  respects,  the  apparatus  and  geometry  were  identical  with  that  of 
Task  1. 

2.2.2  Stationary  Reflections 

Figure  10  is  a  photograph  of  reflections  from  the  stationary  (i.e.,  non¬ 
rotating)  rod.  The  exposure  time  was  1/60  sec.  The  characteristic  appearnace 
of  the  reflections'  diffraction  patterns  results  from  the  hexagonal  shape  of 
lead  carbonate  platelets. 

The  relatively  poor  contrast  in  the  photo  is  due  to  the  strong  background 
scattering  from  the  PMMA  rod.  This  problem  becomes  even  more  severe  with  a 


Figure  10.  Photograph  of  reflections  from  lead  carbonate  mirrors 
embedded  in  the  (stationary)  PMMA  rod. 
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rotating  rod,  where  the  short-lived  reflections  sweep  past  a  nearly  constant 
background,  decreasing  the  signal  to  background  ratio  still  further.  An  effec¬ 
tive  solution  to  this  problem  was  to  reduce  the  duty  cycle  of  the  chopper  by 
covering  5/6  of  the  chopping  wheel  apertures.  This  reduced  the  number  of 
reflection  spots  without  affecting  the  brightness  of  those  that  remain  while 
reducing  the  total  background. 

2.2.3  Task  2  Results  -  The  Rotating  Rod 

The  rod  was  set  in  motion  at  rotation  rate  of  8  Hz.  The  chopping  rate  was 
800  Hz  and  several  photographs  were  made.  The  photograph  included  here  as 
Figure  11  is  representative  of  the  data  obtained.  Several  valid  trajectories, 
defined  as  those  which  have  minimum  of  four  points,  were  identified  and 
labeled  as  A-M  from  film  bottom  to  top.  The  unique  shapes  of  the  reflections 
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aid  in  trajectory  identification  by  eye  and  generally  a  fixed  feature  was  sel¬ 
ected  to  aid  digitization  (centroid,  edge,  etc.)-  Vorticity  values  calculated 
from  these  trajectories  are  compared  with  expected  results  in  Table  3. 


TABLE  3. 


Vorticity  components  calculated  from  photo  trajectories  A-M. 


Trajectory  (No. 

of  Points) 

“z 

“x 

“y 

Bottom  of  Film 

A(5) 

-101.4 

22.9 

-19.7 

B(5) 

-  94.1 

5.7 

1.8 

C(6) 

-104.3 

37.0 

-26.9 

D(6) 

-  94.6 

10.0 

-  3.3 

E(8) 

-  99.5 

23.2 

-20.0 

F(4) 

-  99.5 

38.6 

-32.2 

G(5) 

-  96.2 

17.9 

-15.0 

H(4) 

-100.6 

-  5.8 

1(7) 

-103.1 

19.7 

J(6) 

-  96.8 

20.2 

K(6) 

-  95.4 

18.3 

-12.0 

L(8) 

-  99.6 

10.5 

-  5.6 

Top  of  Film 

M(5) 

-100.3 

16.8 

-10.2 

MEAN 

-  98.9 

18.1 

-11.9 

ST.  DEV. 

3.1 

11.4 

11.0 

OR:  8.1  Hz,  9S  =  12.4°,  <f>s  =  146.7° 


EXPECTED  RESULT:  8.0  Hz,  0S  =  12.0°,  <f>s  =  140.0° 
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2.2.4  Task  2  Discussion 


Vorticity  measurements  tor  many  VOP  mirrors  distributed  throughout  a  finite 
volume  have  been  obtained  simultaneously  in  a  single  exposure.  However,  sev¬ 
eral  features  of  this  data  need  to  be  carefully  examined.  Since  the  number  of 
points  in  each  trajectory  of  Figure  11  are  comparable  to  those  of  Figure  7  of 
Task  1,  and  the  vorticity  components  are  also  comparable,  we  anticipate  the 
standard  deviations  of  the  vorticity  components  in  Table  3  should  be  similar 
to  the  uncertainties  in  vorticity  components  in  Table  1.  That  is,  the  inher¬ 
ent  precision  of  the  vorticity  measurement  in  Task  2  should  be  the  same  in 
Task  1.  Apparently  it  is  not.  Goodness  of  fit  statistics  in  Task  2  were  in 
the  range  of  1  <  X^v  ^  20,  indicating  sources  of  error  other  than  from  digiti¬ 
zation.  The  reason  for  this  discrepancy  is  that  there  is  more  information  in 
these  data  than  has  been  accounted  for  in  our  analysis.  As  described  below, 
some  of  this  information  is  interpreted  as  measurement  error,  which  could  be 
eliminated  by  more  sophisticated  analysis. 

These  "errors"  have  three  distinct  sources.  Digitization  error  has  been 
discussed  in  Section  2.1.5  and  is  a  fundamental  limitation.  The  remaining 
effects  are  due  to  refraction,  particle  position,  and  particle  velocity. 
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Refractive  Effects.  The  film  was  not  in  intimate  contact  with  the  Plexiglas 
block  as  suggested  by  Figure  9.  There  was,  in  fact,  a  4  mm  air  gap  resulting 
in  refraction  at  the  air/plastic  interface  causing  the  effect  indicated  in 
Figure  12.  Rays  which  passed  through  the  center  of  the  block  (normal 
incidence)  were  projected  onto  the  film  plane  correctly.  However,  rays  which 
were  incident  on  the  Plexiglas  surface  at  a  non-normal  angle  were  refracted  to 
strike  the  film  plane  at  an  even  larger  non-normal  angle.  The  greater  the 
angle  of  incidence,  the  greater  the  spot  displacement.  This  introduced 
curvature  artifacts  into  the  trajectories,  which  become  more  severe  near  the 
film  perimeter.  A  straightforward  correction  for  this  effect  was  included  in 
the  data  reduction  program. 
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REFRACTION  DISTORTED 
TRAJECTORIES 

TRUE  TRAJECTORIES 

A-8281 

Figure  12.  Effects  of  refraction  in  air  gap 

(greatly  exaggerated  in  the  figure). 

In  addition,  due  to  a  slight  refractive  index  mismatch  between  DBT  and  our 
PMMA  rod,  there  was  additional  refraction  at  the  PMMA  rod/DBT  surface.  Though 
the  angular  deviations  from  this  source  were  very  small,  the  relatively  long 
path  between  the  rod  and  the  film  resulted  in  spot  displacements  sufficiently 
large  to  affect  trajectory  shape.  Also,  because  the  interface  is  cylindrical 
and  the  origin  of  the  reflections  within  unknown,  this  effect  is  essentially 
uncorrectable.  It  is,  however,  an  artifact  of  this  particular  experiment  and 
we  emphasize  that  it  would  not  exist  under  actual  flow  conditions  with  actual 
VOP  particles. 


Particle  Position  and  Velocity  Effects.  The  trajectory  analysis  used  through¬ 
out  this  work  assumes  that  all  reflections  emanate  from  the  origin  of  the 
coordinate  system  in  Figure  2.  However,  due  to  the  non-zero  diameter  of  the 
rotating  rod,  reflections  actually  originate  offset  from  the  origin  and  they 
translate  as  the  rod  rotates.  These  features  result  in  apparent  errors  in 
vorticity  determinations.  Knowledge  of  the  VOP  mirror's  position  for  each 
spot  on  the  trajectory  would  re-establish  the  inherent  precision  of  vorticity 
measurement  seen  in  Task  1,  with  the  additional  benefit  of  providing  the 
particle  velocity.  Conversely,  it  could  be  said  that  to  obtain  the  vorticity 
measurement  accuracy  of  Task  1,  we  must  either  measure  the  particle  velocity 
or  restrict  the  size  of  the  measurement  volume.  However,  all  of  the  infor¬ 
mation  needed  to  deduce  the  particle  velocity  can  be  encoded  in  the  reflection 
trajectories;  an  independent  measure  of  velocity  is  not  required.  Demonstra¬ 
ting  this  capability  was  the  goal  of  Task  3. 
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2.3  TASK  3  -  LOCATING  REFLECTION  ORIGINS  (one-dimensionai) 

To  study  the  shape  and  evolution  of  a  coherent  structure  by  measuring  its  char¬ 
acteristic  vorticity  field,  the  precise  position  and  time  from  within  a  large 
volume  at  which  a  local  vorticity  vector  is  measured  must  be  known.  In 
addition,  we  learned  in  Task  2  that  reflection  origins  must  be  known  to 
achieve  maximum  precision  in  vorticity  measurement.  For  these  reasons,  Task  3 
was  designed  to  demonstrate  the  feasibility  of  encoding  reflection  origin 
information  in  the  reflected  light  itself. 
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2.3.1  Task  3  approach 
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If  the  sampled  volume  is  illuminated  with  two  coincident  beams  of  different 
color,  the  reflections  will  of  course  contain  the  two  colors  as  well.  Then, 
if  known  spatial  intensity  gradients  are  imposed  on  one  or  both  of  the  colors 
along  the  z-direction  for  example  as  in  Figure  13,  it  is  possible  to  determine 
the  Z  coordinate  by  observing  the  ratio  of  reflected  intensity  of  the  two 
colors. 


For  intensity  gradients  in  color  1  and  color  2  respectively  of: 


II  =  f(z)  (e.g.  az  +  b),  l£  =  g(z), 


(12) 


v. 


.7, 


\ 


3 


the  ratio  of  reflected  intensity  at  position  Z  is: 

l!/I2  =  f(z)/g(z)  s  h(z) ,  (13) 


Figure  13.  Coincident  beams  of  different  color  or  polarization. 
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where  h(z)  is  a  well-determined  monotonic  function  of  Z.  The  measured  I\/l2 
ratio  corresponds  to  a  unique  Z  value.  This  procedure  can,  in  principle,  be 
extended  to  multiple  dimensions  but,  as  discussed  later,  other  methods  in 
conjunction  with  the  one-dimensional  measurement  may  eliminate  the  need  for 
more  than  two  colors. 


We  had  proposed  to  perform  this  measurement  by  separating  the  5145A  and  4480A 
lines  emitted  by  our  Argon-Ion  laser,  imposing  an  intensity  gradient  on  one 
color,  and  then  recombining  the  two  beams  prior  to  illuminating  the  PMMA  rod. 
The  reflections  were  to  be  captured  on  color  photographic  film,  which  would 
have  been  analyzed  by  microdensitometry  to  yield  the  color  intensity  ratios. 
Unfortunately  we  had  not  realized  at  the  time  that  these  two  colors  are  not 
sufficiently  distant  from  one  another  in  wavelength  to  be  easily  discriminated 
in  color  photographic  film.  Since  no  other  color  or  laser  was  immediately 
available,  we  chose  an  alternate  approach  which  demonstrated  the  validity  of 
the  method,  albeit  in  a  somewhat  less  colorful  manner,  when  used  in  conjunc¬ 
tion  with  modern  electro-optic  detectors  and  processors. 


2.3.2  Apparatus 


The  simple  demonstration  of  this  measurement  employed  a  single  color  of  light 
and  two  successive  stationary  reflections,  illuminated  by  beams  having  two 
different  imposed  intensity  gradients.  A  digital  video  system  was  used  to 
capture  and  compare  images  of  the  reflections.  Figure  14  shows  the  Task  3 
optical  configuration. 


The  incident  beam  was  colinear  with  a  pivot  about  which  the  optical  apparatus 
(but  not  the  rod  or  the  chamber)  turns.  This  allowed  adjustment  (via  the 
micrometer  screw)  of  beam  position  along  the  rod's  axis  without  altering  the 
optical  path.  The  beam  was  expanded  and  passed  through  a  rectangular 
aperture,  thereby  providing  a  nearly  uniform  intensity  "top-hat"  beam  with 
sharp  edges.  The  top-hat  beam  passed  through  a  Newport  Corporation  linear 
gradient  neutral  density  filter,  which  imposed  a  roughly  linear  intensity 
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Figure  14. 


Task  3  optical  configuration  for  intensity  gradient 
position  measurement. 


gradient  from  top  to  bottom  (the  z-axis,  normal  to  page  in  the  figure).  This 
gradient  filter  could  be  inverted  to  provide  an  intensity  gradient  in  the 
opposite  sense.  The  beam  was  then  condensed  to  desired  cross-section  and 
directed  through  the  DBT  chamber  to  the  mirror-bearing  rod.  The  reflections 
on  the  Plexiglas  diffuser  plate  were  then  viewed  with  a  Sierra  Scientific  CCD 
camera  equipped  with  a  50  mm  lens.  The  camera  output  was  coupled  to  a  Data 
Translation  Model  DT2851,  Frame  Grabber  and  Ceii  aq  286  computer,  which  enabled 
capture  and  digitization  of  individual  camera  frames.  Each  frame  comprised 
512  x  512  individual  pixels  digitized  to  8-bit  accuracy.  The  amplitude  of  the 
digitized  signal  at  each  pixel  is  linearly  proportional  to  the  energy  incident 
on  that  pixel  during  the  exposure  time,  plus  a  fixed  offset. 


2.3.3  Experimental  Procedure 
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Direct  Position  Measurements.  In  order  to  test  the  intensity  gradient  tech¬ 
nique,  the  Z  positions  of  individual  mirrors  first  were  measured  directly  and 
then  compared  to  the  positions  predicted  by  the  gradient  method.  The  direct 
measurements  were  accomplished  by  using  the  micrometer  screw  to  raise  and 
lower  the  rectangular  beam  along  the  mirror  bearing  rod's  axis  (i.e.,  the  z- 
axis).  From  a  standard  starting  position,  the  beam  was  lowered  until  a 
particular  reflection  vanished  from  the  diffuser  plate,  then  raised  from  this 
position  until  it  vanished  once  more.  This  gave  the  beam  height.  The  beam 
was  then  lowered  to  the  standard  position  leaving  the  particle  at  a  known  Z 
position  within  the  incident  beam.  In  this  way  many  individual  reflections 
were  measured  and  catalogued  as  the  motor  shaft  was  burned  by  hand. 

Intensity  Gradient  Measurements.  To  compare  the  direct  position  measurements 
with  the  intensity  gradient  measurement,  two  images  were  captured  for  each 
reflection  or  group  of  reflections  within  the  CCD  camera  field  of  view.  The 
first  image  had  the  approximately  linear  intensity  gradient  oriented  so  that 
the  brightest  part  of  the  incident  beam  was  near  the  top.  This  is  abbreviated 
as  the  gradient  "up"  image.  The  gradient  filter  was  then  inverted  to  capture 
the  second  gradient  "down"  image. 

To  analyze  these  images,  subtraction  of  the  fixed  background  is  required,  as 
is  a  calibration  of  the  intensity  gradient.  To  measure  the  background  signal, 
the  rod  was  rotated  to  a  position  where  no  reflections  could  be  seen  in  the 
field  of  view  and  only  the  bright  background  was  visible.  This  scene  was 
recorded  gradient  "up"  and  "down".  The  calibration  was  accomplished  by  a  simi¬ 
lar  method.  The  DBT  chamber  and  rod  were  removed,  and  the  incident  beam  was 
projected  directly  onto  the  diffuser  screen.  Images  of  it  were  also  recorded 
gradient  "up"  and  "down". 

Image  analysis.  A  video  processing  routine  was  written  to  perform  background 
subtraction  and  divide  "up"  frames  by  "down"  frames.  The  resulting  ratios 
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were  found  to  be  quite  uniform  over  an  individual  reflection's  image. 

Figure  15  shows  such  a  processed  image.  The  frame  division  has  been  performed 
only  over  the  rectangular  area  on  the  spot  (at  upper  left)  since  the  signal/ 
noise  ratio  degrades  near  the  reflection's  edge.  An  intensity  profile  across 
the  reflection  (along  the  line  visible  inthe  photo)  is  also  shown.  The  values 
are  one-hundred  times  the  actual  intensity  ratio. 

Similar  processing  of  the  calibration  image  provided  a  measure  of  intensity 
ratios  at  every  position  in  the  incident  beam.  The  calibration  and  its  corres 
ponding  profile  are  shown  in  Figure  16.  To  establish  the  z-coordinate  of  a 
particular  reflection  origin,  it  remains  only  to  find  the  position  along  the 
calibration  curve  of  Figure  16  where  the  intensity  ratio  is  the  same  as  that 
observed  in  the  reflection. 


Figure  15.  False  color  enhanced  video  display  of  gradient  "up"  frame 

(background  subtracted)  after  division  by  the  corresponding 
"down"  frame  within  the  rectangular  area  on  the  reflection  at 
upper  left. 


Figura  16.  Calibration  band.  Division  has  been  performed  over  the 
directly  imaged  incident  beam  cross-section. 

2.3.4  Task  3  Results 

Figure  17  is  a  plot  of  the  direct  measurements  against  the  intensity  gradient 
measurements.  The  calibration  curve  was  also  employed  to  convert  uncertainty 
in  reflection  ratios  to  uncertainty  in  particle  positions.  Position  errors 
result  from  small  variations  in  gradient  uniformity  or  differences  in  the 
optical  path  between  successive  images.  The  intensity  gradient  method  gives 
~200  urn  precision  over  the  3.5  mm  beam.  This  is  less  than  10  VOP  particle 
diameters . 


2 . 4  CONCLUSIONS 


In  this  program,  experiments  were  performed  that  demonstrated  the  feasibility 
of  expanding  the  extant  vorticity  optical  probe  technology  to  include  three 
component  vorticity  measurement  and  volumetric  applications.  Precise  three 
component  vorticity  measurement  has  been  demonstrated  at  a  single  point  and  at 
many  points  within  a  volume.  However,  it  has  been  found  that  the  best 
accuracy  is  achieved  when  the  location  of  the  reflection  origin  is  well  known 
for  each  point  (spot)  on  the  trajectory.  A  possible  method  for  determining 
these  locations  has  been  demonstrated,  using  state-of-the-art  digital  video 
technology. 

Knowledge  of  particle  positions  along  a  reflection  trajectory  (i.e.  as  a 
function  of  time)  is,  of  course,  equivalent  to  particle  velocity  measurement. 
Though  the  feasibility  of  measuring  one  position  coordinate  has  been  demon¬ 
strated,  its  extension  to  more  than  one  dimension  may  not  be  straightforward. 
The  difficulty  is  in  obtaining  multiple  colors  of  sufficient  intensity  to 
enable  detection  of  reflections  rotating  at  speeds  as  much  as  100  times  those 
used  in  this  Phase  I  program.  It  is  conceivable  that  this  information  can  be 
obtained  from  the  reflection  trajectories  by  performing  a  more  complete  mathe¬ 
matical  analysis.  Specifically,  the  analysis  of  Section  1.2.1  could  be  modi¬ 
fied  to  include  the  six  additional  unknowns  X0,  Y0,  Z0,  U,  V,  W.  However, 
preliminary  examination  of  the  full  (nine  unknowns  with  w)  description  of 
trajectories  indicates  that  many  of  these  variables  may  be  too  strongly 
coupled  to  be  determined  from  a  single  trajectory.  Furthermore,  they  do  not 
enter  in  the  simple  linear  fashion  that  so  greatly  simplified  the  vorticity 
analysis.  It  would  be  preferable  to  design  experiments  which  provide  the 
auxiliary  data  to  make  such  an  analysis  possible.  A  potential  approach,  shown 
in  Figure  18,  employs  coincidence  detection  on  two  opposed  detectors. 
Oppositely  directed  beams  produce  spot  pairs  on  the  detectors;  the  reflection 
origin  must  lie  on  the  line  joining  these  spots.  An  independent  depth  deter¬ 
mination  (along  y  in  the  figure)  using  a  version  of  the  intensity  gradient 
method,  for  example,  would  give  complete  position  and  velocity  information. 
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Figure  18.  Coincidence  method  for  complete  vorticity  and  velocity 
determination. 


Also  the  similarity  of  the  patterns  on  the  detectors  (with  simple  corrections 
for  the  mean  velocity)  would  facilitate  automated  trajectory  identification 
and  data  reduction. 


For  non-volumetric  applications,  the  VOP  with  three  component  capability  is 
itself  a  valuable  instrument  which  may  have  commercial  appeal.  The  principal 
obstacle  to  its  wide-spread  applicability  remains  the  lack  of  water  compati¬ 
bility.  We  have  studied  this  problem  and  believe  that  there  are  several 
possible  solutions.  The  most  appealing  is  to  manufacture  VOP  particle  from 
one  of  a  variety  of  low  refractive-index  fluorinated  oils  or  plastics.  Upon 
solving  this  problem,  we  believe  that  a  prototype  (single  point)  VOP  is 
quickly  realizable. 

We  hope  to  have  the  opportunity  to  address  and  overcome  these  issues  in  a 
Phase  II  program.  The  first  objective  of  Phase  II  would  be  to  manufacture  VOP 
particles  which  can  be  used  in  common  laboratory  water  flow  tunnels.  Upon 
achieving  this  goal,  we  would  proceed  with  development  of  a  prototype  commer¬ 
cial  system  able  to  measure  the  vorticity  vector  at  a  single  point  in  the 
flow.  A  small  water  flow  system  would  be  assembled  for  use  in  testing  the 
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prototype.  Concurrently,  we  would  continue  the  development  of  a  volumetric 
measurement  capability,  addressing  the  details  of  three-coordinate  position 
measurement  and  of  tracking  multiple  trajectories  from  rapidly  rotating  mirror¬ 
bearing  spheres. 
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An  optical  method  for  the  direct  measurement  of  vorticity  in  liquid  flows  is  described. 
At  the  present  state  of  development  it  is  capable  of  responding  to  vorticity  fluctuations 
with  a  response  time  of  about  1  msec  and  a  spatial  resolution  of  better  than  50 /an. 
Small  spherical  particles  suspended  in  the  flow  rotate  with  angular  velocity  accurately 
equal  to  half  the  local  vorticity;  thus  measurements  of  the  rotation  rates  of  such 
particles  indicate  the  vorticity.  Transparent  spherical  particles  of  less  than  50  /tm 
diameter,  each  containing  embedded  planar  crystal  mirrors,  have  been  developed  for 
this  purpose  and  are  suspended  in  a  refractive-index-matched  liquid.  Measurements  of 
the  times  required  for  laser  reflections  from  the  mirrors  to  rotate  through  the  small  angle 
defined  by  a  pair  of  slits  yields  the  rotation  rate,  and  thus  the  vorticity.  Production  and 
physical  properties  of  the  probe  particles  are  reported.  Theoretical  capabilities  and 
limitations  of  the  method,  including  accuracy,  spatial  and  temporal  resolution,  data 
rate,  and  background  noise  are  calculated  and  found  to  be  coupled  to  the  optical 
geometry  and  flow  field.  Analysis  yields  procedures  for  selective  optimization  of  each 
parameter  as  dictated  by  the  particular  application.  Measurements  of  steady -state, 
laminar,  two-dimensional  Poiseuille  flows  demonstrate  the  effectiveness  of  the 
technique  and  confirm  theoretical  predictions. 


1.  Introduction 

Although  the  vorticity  field.  u>(x,t)  =  V  x  v(x,f),isanessentialproperty  of rotational 
and  turbulent  flow,  requisite  measurements  of  the  vorticity  transcend  experimental 
capability.  Modem  techniques  of  hot-wire  andhot-filmanemometryand  laser-DoppIer 
anemometry  can  measure  components  of  the  velocity  field  v(x,  t)  at  a  limited  number 
of  points.  The  velocity  components  are  ordinarily  analysed  by  finite  differences  or 
application  of  Taylor’s  frozen-flow  hypothesis  to  calculate  the  vorticity  (see,  for 
example,  Eckelmann  e<  al.  1977;  Foss  1977;  Willmarth  &  Bogar  1977),  but  error 
magnification  by  the  necessary  differentiations  limits  the  accuracy  and  spectral  range. 
Spatial  resolution  of  the  method  has  usually  been  limited  to  about  0- 1  cm  by  hot-wire 
probe  geometry  although  Willmarth  &  Bogar  (1977)  have  recently  reported  a 
spatial  resolution  of  100 /tm.  In  principle  flow  perturbations  can  be  avoided  and  the 
spatial  resolution  can  be  slightly  improved  by  laser-Doppier  anemometry  but  the 
complexity  of  these  simultaneous  optical  measurements  has  apparently  deterred  its 
application  to  vorticity  measurement. 

The  ideal  system  for  vorticity  measurement  should  meet  these  criteria:  Spatial 
resolution  exceeding  the  inner  scale  of  turbulence  ij  (rj  ~  50 /nn  in  common  liquid 
flows),  temporal  resolution  exceeding  the  shortest  expected  vorticity  fluctuation  time. 
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and  data  acquisition  rates,  spectral  range  and  measurement  accuracy  sufficient  for 
calculation  of  vorticity  correlation  functions  and  higher  moments  of  vorticity. 

The  potential  of  such  measurements  of  vorticity  has  motivated  a  disparate  approach 
to  this  difficult  measurement  problem.  Encouraged  by  the  successful  development  of 
a  method  of  optical  scattering  from  anisotropic,  submicroscopic  particles  for  direct 
measurement  of  the  strain-rate  tensor  with  excellent  spatial  resolution  (Johnson  & 
Webb  1972),  a  new  optical  method  for  direct  measurement  of  vorticity  distributions  in 
liquids  with  50 /im  spatial  resolution  has  been  developed  and  is  described  here.  It  is 
based  on  the  propensity  of  small  spherical  particles  suspended  in  flowing  fluid  to 
accurately  track  the  local  vorticity  u>  by  rotation  with  angular  velocity  SI  -  Jw. 
Reflections  from  a  dispersion  of  small  transparent  plastic  spheres  each  containing  one 
or  more  flat,  highly  reflecting  mirrors  provide  the  desired  indicator  of  local  vorticity. 
To  measure  vorticity,  a  few  volume  per  cent  of  these  particles  are  added  to  a  trans¬ 
parent  liquid  of  matching  refractive  index.  A  laser  beam  illuminating  a  small  volume 
of  the  flow  to  be  sampled  is  reflected  by  the  rotating  particles.  The  rotational  velocity 
of  each  reflected  beam,  which  equals  the  local  vorticity  at  the  reflecting  particle,  is 
measured  by  observing  the  transit  time  of  the  reflected  beam  across  a  small  angle 
defined  by  a  pair  of  closely  spaced  slits.  Each  such  transit  provides  a  datum  for  prompt 
on-line  calculation  of  the  instantaneous  local  vorticity  that  is  stored  for  analysis  of  the 
vorticity  distribution.  This  measurement  system  is  here  called  the  vorticity  optical 
probe  or  VOP  for  brevity. 

The  vorticity  measurement  system  is  described  in  §2  of  this  paper.  There  the 
processs  of  manufacture  and  the  physical  properties  of  the  probe  particles  are  reported 
in  some  detail.  This  is  followed  by  a  description  of  the  data  collection  system  designed 
for  application  to  laboratory  scale  wakes,  mixing  layers  and  boundary  layers  in 
channels  at  modest  Reynolds  numbers.  The  characteristics  of  the  detected  reflexions, 
which  tend  to  limit  the  precision  of  the  technique,- depend  on  both  reflecting  particle 
properties  and  optical-system  geometry.  The  optical  geometry  its  and  effect  on 
measurement  performance  are  discussed  and  the  signal  acquisition  and  data-processing 
procedure  are  described. 

Section  3  presents  a  theoretical  discussion  of  basic  experimental  capabilities  and 
limitations  of  the  method  including  accuracy,  spatial  and  temporal  resolution,  data 
rate  and  background  noise,  and  their  dependence  on  experimental  parameters.  To 
describe  how  the  interplay  of  inhomogeneity  of  flow  with  the  optical  geometry  and 
the  optical  imperfections  limits  the  available  vorticity  measurement  resolution  a 
system-resolution  function  is  defined. 

Section  4  describes  experimental  tests  of  the  measurement  system  by  application  to 
steady-state  Poiseuille  flows.  Measured  vorticity  distributions  confirm  theoretical 
expectations  and  demonstrate  the  effectiveness  of  the  vorticity  optical  probe. 

Finally  in  §  5  the  properties  of  the  vorticity  optical  probe  system  are  summarized 
and  prospects  for  applications  and  future  development  are  discussed. 


& 

< 

) 

I 

§ 

S 


i 


£ 


I 

*  J 

# 

* 

f 

*.j« 


;•■ 


2? 


I 


§ 


7T 


T7T"TT* 


P!P*','J 


I 


Direct  measurement  of  vorticity  by  optical  probe  175 

2.  Vorticity  measurement  system 

2.1.  Fundamental  concepts 

The  vorticity  measurement  system  comprises:  (1)  the  spherical  probe  particles  con¬ 
taining  planar  mirrors;  (2)  refractive- index- matched  fluid;  (3)  optical  system;  (4)  data 
acquisition  and  analysis  system. 

Each  aspect  of  the  system  will  be  discussed  in  turn.  First,  however,  the  fundamental 
hypothesis  that  small  spherical  particles  follow  the  flow  must  be  established. 

Chwang  &  Wu  (1974)  have  calculated  the  couple  due  to  local  vorticity  acting  on 
small  particles  of  various  shapes  with  results  that  agree  with  the  early  calculation  of 
Jeffrey  (1922)  for  spheroids.  Application  of  their  results  shows  that  the  spherical  probe 
particles  faithfully  follow  flow  fluctuations  that  occur  on  a  time  scale  slower  than 
about  1  ms.  The  calculations  are  valid  when  the  Reynolds  number  based  on  particle 
radius  r  and  velocity  u  relative  to  the  surrounding  fluid  is  much  less  than  unity. 
Stokes’  law  suggests  that  relative  velocity  will  relax  exponentially  with  relaxation 
time 

Tu  =  ik*p//*.  (2.1) 

where  r  is  the  sphere  radius,  p  the  particle  density,  and  p  the  dynamic  viscosity.  With 
appropriate  data  for  the  probe  particles,  p  =  l-2gcm-3,  p  -  0-005 g cm-1  s-1  and 
r  =  25  x  10~4  cm,  the  relative  velocity  relaxation  time  is  ru  ~  0-4  ms.  Hence,  the 
calculation  of  Chwang  &  Wu  is  valid  for  the  probe  particles  for  time  scales  >  1  ms. 

Chwang  &  Wu’s  equation  (11)  gives  the  torque  on  a  sphere  of  radius  r  rotating  at 
angular  velocity  Q  relative  to  the  surrounding  fluid  as  M  =  —  8 npr3^!.  Thus  the 
relative  rotational  relaxation  time  is 

Tr  =  T  h^p/p.  (2.2) 


For  the  probe  particles  r,  ~  0-lms. 

The  shortest  time  scale  in  turbulent  flow  is  the  Kolmogorov  micro-scale  of  time 
(Tennekes  &  Lumley  1972) 

Tk  ((/w)i?i  (2.3) 

where  Rt  is  the  turbulence  Reynolds  number,  l  the  integral  length  scale,  and  u  the 
turbulent  velocity.  Extreme  values  in  laboratory  water  tunnel  flows  are  l  ~  0-3  cm, 
u  ~  7cms-1,  and  R ,  ~  200  (Frenkiel,  Klebanoff  &  Huang  1979),  giving  rk  —  3ms. 
Therefore  the  probe  particles  should  readily  follow  rotational  and  translational 
fluctuations  in  anticipated  turbulent  liquid  flows. 

This  line  of  analysis  also  provides,  following  Jeffrey  (1922),  an  estimate  of  the 
possible  measurement  error  due  to  non-spherical  particles  that  might  be  formed  by 
damage  or  adhesion  of  spherical  probe  particles.  Non-spherical  particles  do  not  rotate 
with  the  same  angular  velocity  as  the  fluid.  Their  long  axes  tends  to  align  with  the 
principal  strain  rate  in  anon-uniform  flow.  Approximating  an  adherent  pair  of  spheres 
by  an  ellipsoid  of  major  axis  4r  yields  a  variable  rotation  rate  of  the  form 
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Figure  1.  Photomicrograph  of  basic  lead  carbonate  crystals  in  reflected  light. 


where  dv/dy  is  the  maximum  magnitude  of  the  local  velocity  gradient.  The  rotation 
rate  Qp  of  this  ellipsoid  varies  from  the  rate  Q  for  spherical  particles  within  the  limited 
range 

<  Qp  <  fil.  (2.5) 

Thus  defective  probe  particles  would  contribute  a  broadened  background  to  vorticity 
measurements. 

Rotation  of  the  probe  spheres  should  not  be  perturbed  by  inter-particle  interaction 
and  they  should  not  perturb  the  fluid  flow.  Jeffrey  (1922)  showed  that  the  effective 
viscosity  of  the  bulk  fluid  will  be  increased  by  a  factor  of  1  +  2-5  V  by  the  addition  of 
a  volume  fraction  V  of  solid  spherical  particles.  Thus  the  probe  particle  number  density 
n  should  be  limited  so  that  the  volume  fraction  of  particles  is  small,  that  is  f  ni^n  <g  1 . 
For  spheres  of  25 /tm  radius  this  requires  n  1-5  x  107  cm-3.  A  more-than-ample  con¬ 
centration,  10*  particles  cm"1,  increases  the  viscosity  by  only  10  %. 


2.2.  Probe  particles 

The  preceding  section  indicates  that  spherical  particles  of  radius  r  <  25  fim  and 
moderate  density  accurately  sample  the  local  vorticity.  They  rotate  at  an  angular 
velocity  Q  equal  to  half  the  local  vorticity  w  (Batchelor  1967,  §§2.3,  2.4).  Uniform 
dispersion  of  the  particles  requires  that  they  not  stick  together  and  that  the  particle 
density  approximate  that  of  the  fluid.  To  minimize  light  scattering  the  refractive  index 
of  the  particles  must  closely  match  that  of  the  fluid.  The  essential  development 
enabling  this  optical  method  for  vorticity  measurement  is  the  process  for  production 
of  large  quantities  of  transparent  spherical  particles  containing  planar  mirrors.  The 
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Figure  3.  Histogram  of  particle  diameter  distribution. 

reflexion  of  a  laser  beam  from  a  mirror-bearing  sphere  rotates  with  angular  velocity 
equal  to  the  local  vorticity  o». 

The  following  subsections  describe  the  manufacture  and  characteristics  of  the  probe 
particles. 

(a)  Particle  manufacture 

Methylmethacrylate  monomer,  a  volatile,  low- viscosity  liquid,  is  easily  ‘mass’ 
polymerized  (Rodriguez  1970.  Redfam  &  Bedford  1960)  by  adding  a  small  amount 
( 0- 1  %  by  weight)  of  initiator  (usually  benzoyl  peroxide)  and  heating  for  several  hours 
A  clear,  hard,  plastic  solid  (commercially  known  as  Lucite  or  Plexiglas)  results.  If  the 
non- water-soluble  monomer  is  suspended  in  a  continuously  stirred  aqueous  solution, 
it  breaks  up  into  small  spherical  ‘  pearls  ’  and  the  same  chemical  polymerization  occurs. 
This  procedure  is  known  as  suspension  or  pearl  polymerization  since  the  pearls  retain 
their  shape  when  hardened. 

The  suspension  fluid  consists  of2g/l  polyvinyl  alcohol  (PVA),  l-5g/leacb  of  sodium 
chloride  (NaCl)  and  ammonium  thiocyanate  (NH4SCN)  in  doubly  distilled  water.  PVA 
prevents  the  liquid  pearls  from  coalescing  once  they  are  separated,  NaCl  helps  to 
prevent  partially  polymerized  pearls  from  sticking  together,  and  NH4SCN  inhibits  the 
formation  of  emulsion  polymer  (Rodriguez  1970;  Hopff,  Liissi  &  Gcrspacher  1964). 
Approximately  150  ml  of  monomer  plus  initiator  is  added  to  850  ml  of  solution,  and 
then  whipped  for  5  s  in  a  Montgomery  Ward  pulse  action  blender.  The  suspension  is 
then  poured  into  a  one-litre  volumetric  flask,  heated  to  65-70  °C  and  magnetically 
stirred  on  a  Corning  heater /stirrer  for  about  five  hours.  The  spheres  harden  and, 
when,  dried  form  about  100  ml  of  slippery  powder. 

Lead  carbonate  ‘mirrors’  are  easily  incorporated  in  the  particles  before  they 
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polymerize  by  adding  7-5 ml  of  Mearlmaid  Nacromer  ZTX-B  (supplied  by  Mearl  Co.) 
to  the  monomer  before  whipping.  The  mirrors  are  hexagonal  platelets  with  widths  of 
8-30 pm,  thickness  of  0-01  pm  and  refractive  index  of  2-2  (see  figure  1).  The  mirrors 
generally  remain  fully  enclosed  by  the  plastic  spheres  throughout  the  process  and  do 
not  interfere  with  the  polymerization  chemistry. 

(b)  Particle  characteristics 

The  resulting  spheres  have  a  densityofabout  l-2gcm-3andarefractive  index  of  1-49. 
Microscopy  shows  a  diameter  range  of  about  5-200 pm.  The  smallest  spheres  contain 
no  mirrors,  the  larger  ones  contain  several.  Since  only  spheres  smaller  than  a  Kolmo¬ 
gorov  microscale  (Tennekes  &  Lumley  1972)  of  roughly  50pm  that  contain  mirrors  are 
desired,  gravitational  sedimentation  (Cadle  1955)  is  used  to  narrow  the  size  range. 
Figures  2  (a,b)  show  phase-contrast  photomicrographs  of  finished  particles  imbedded 
in  an  index-matching  fluid  and  illuminated  by  transmitted  and  reflected  light  respec¬ 
tively.  Measurements  of  many  similar  pictures  determined  the  final  particle  diameter 
distribution  shown  in  figure  3.  The  average  diameter  is  21  pm.  It  appears  that  there 
are  on  the  average  about  two  mirrors  per  sphere. 

To  check  for  flocculation  and  non-spherical  or  otherwise  defective  particles  a  micro¬ 
scopic  flow  system  was  constructed  in  which  the  particles  could  be  observed  under 
actual  flow  conditions.  It  was  found  that  a  small  minimum  shear  rate,  which  is  easily 
surpassed  in  any  laboratory  flow  system,  sufficed  to  keep  the  spheres  dispersed.  The 
observations  also  show  that  approximately  3  %  of  the  particles  have  defects  such  as 
fractures  or  non-spherical  form.  An  additional  3-4  %  appear  to  be  stuck  together  in 
inseparable  pairs  or  occasional  triplets. 

2.3.  Fluid  selection 

The  particles  must  be  suspended  in  a  fluid  having  the  following  properties :  ( 1 )  refractive 
index  which  nearly  matches  that  of  the  spheres  (1.49)  to  prevent  scattering  from  their 
surfaces  that  would  limit  the  optical  depth;  (2)  no  dissolution  or  chemical  reactivity 
with  the  particles  or  components  of  the  flow  system ;  ( 3)  density  near  that  of  the  particles 
to  prevent  rapid  settling  out;  (4)  low  viscosity  to  facilitate  experiments  with  high 
Reynolds  number;  (5)  non-foaming  surface  properties;  (6)  low  vapour  pressure, 
toxicity  and  flammability  for  safety  and  convenience;  and  (7)  low  cost  to  permit 
economical  extension  to  large-scale  apparatus. 

Two  fluids,  dibutyl  phthalate  and  p-cvmene,  were  found  which  satisfy  all  or  most  of 
the  requirements.  The  former  has  a  viscosity  of  0-15gcm_1s-1  which  is  suitable  for 
the  laminar  flow  experiments,  while  the  latter  has  lower  viscosity,  0-005  g  cm-1  s-1, 
but  somewhat  higher  flammability  and  vapour  pressure.  Certain  concentrated  solu¬ 
tions  of  salts  in  water  that  also  satisfy  most  of  these  critera  are  excessively  corrosive. 

2.4.  Detection  system 

The  detection  system  consists  of  two  parts,  the  optical  system  and  the  electronics  for 
data  collection  and  processing.  The  optical  system  is  designed  to  perform  two  func¬ 
tions:  Define  the  sampled  volume  in  the  fluid,  and  define  the  particle  rotation  angle  to 
be  timed  by  the  electronics.  The  sampled  volume  is  defined  by  the  illuminating  laser 
beam  and  a  lens-slit  combination  that  collects  the  reflections  and  limits  the  observed 
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Figure  4.  Optical  detection  system. 

length  of  the  laser  beam  as  will  be  described.  Another  pair  of  slits  with  fixed  separation 
defines  the  rotation  angle.  Light  passing  through  these  slits  is  converted  to  electrical 
pulses  by  photomultipliers  (RCA  7265).  The  pulses  trigger  a  pair  of  threshold  detectors, 
one  of  which  {on)  starts  a  clock  on  a  pulse  from  the  first  PMT  and  the  second  {off)  stops 
the  clock  on  a  pulse  from  the  second  PMT.  The  elapsed  time  is  read  by  a  computer 
which  immediately  calculates  the  vorticity.  The  accuracy  of  the  data  is  determined 
by  the  regularity  of  the  light  beams  passing  through  the  slits,  which  in  turn  is  deter¬ 
mined  by  the  quality  of  the  reflexions  from  the  particles  and  the  quality  of  the  optical 
components.  In  this  section  the  detection  system  is  described  in  detail  to  provide  a 
basis  for  a  determination  of  the  measurement  uncertainty  function. 

(a)  Optical  system 

(i)  Optical  system  geometry.  Figure  4  shows  a  schematic  diagram  of  the  optical  system 
and  its  relationship  to  the  flow  system  The  incident  beam  enters  the  channel  at  a 
45°  angle  through  a  prism  which  is  refractive-index-matched  to  the  channel  wall  The 
reflexions  from  suitably  oriented  particles  emerge  nearly  perpendicular  to  the  incident 
beam  through  a  second  prism  and  are  imaged  without  magnification  onto  slit  3,  as 
indicated  by  the  ray  tracing  in  figure  4.  The  incident  beam  and  slit  3  define  the  sampled 
volume.  Slits  1  and  2  are  fixed  in  position;  along  with  movable  mirrors  1  and  2  they 
define  the  particle  rotation  angle,  2a,  through  which  a  particle  rotating  in  the  direction 
indicated  must  turn  for  its  reflexion  to  pass  from  slit  1  to  slit  2.  Since  the  reflexions 
rotate  twice  as  fast  as  the  particles,  the  slits  define  an  actual  angle  of  4a.  The  maximum 
height  of  these  slits  in  this  apparatus  is  5  cm,  limited  by  the  size  of  the  photomultiplier 
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Figure  5.  Possible  positions  of  reflected  diffraction  spot  relative  to  slits, 
illustrating  its  truncation. 
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Figure  6.  Segment  of  dual -channel  chart  recorder  trace  showing  output  of  photomultipliers  due 
to  reflexions  from  particles  embedded  in  a  rigidly  rotating  rod. 

anodes  placed  directly  behind  them.  Consequently,  only  about  4  %  of  the  total  number 
of  reflexions  are  intercepted  by  the  optics.  A  more  elaborate  apparatus  with  specially 
shaped  slits  and  larger  detectors  could  increase  the  collection  efficiency  substantially. 
All  of  the  optical  components  are  movable,  as  shown  by  the  double-headed  arrows  in 
figure  4,  to  allow  repositioning  of  the  sampled  volume  without  significantly  changing 
the  reflexion  path  length  of  37  cm  or  the  rotation  angle  2 a. 

(ii)  Properties  of  the  reflexions.  Reflexions  of  parallel  illumination  from  perfect 
mirrors,  approximated  by  disks  of  radius  L,  are  diffraction  broadened  to  Airy  disk 
patterns  with  the  central  disk  spreading  at  an  angle  ft  =  1-22A/2L.  For  the  probe 
mirrors  ft  ~  15°.  For  an  incident  beam  intensity  I,  the  reflected  power  is  IttL1.  The 
observed  irradiance  on  a  plane  a  distance  z  from  the  mirror  is 

Ir  =  InL*/[Tr(ftz )«]  =  2-69/L7AV*.  (2.6) 

Thus  the  intensity  of  the  reflexion  varies  with  the  fourth  power  of  mirror  size .  Additional 
variations  can  arise  from  the  following  factors:  ( I )  particles  move  through  an  incident 
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Figure  7.  Distribution  of  elapsed  times  between  threshold  crossings  of  detector  1  followed  RN 
detector  2,  as  determined  from  data  like  that  of  figure  6.  M 
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Figure  8.  Block  diagram  of  data  analysis  system. 

beam  with  a  gaussian  intensity  profile  while  they  rotate,  thus  changing  reflexion" 
intensity;  (2)  as  shown  in  figure  5,  portions  of  the  reflected  beam  may  be  cut  off  by  the"’ 
top  or  bottom  of  a  slit;  (3)  the  intensity  of  the  ideal  reflexions  is  angle  dependent; 
(4)  mirror  imperfections  as  well  as  differences  in  mirror  diameter  vary  the  intensity’ 
between  particles;  and  (5)  optical  screening  by  separate  particles  can  attenuate  th  * 
reflexion. 

To  characterize  the  reflexions  a  dummy  specimen  consisting  of  particles  imbeddet?^ 
in  a  rigid  rod  of  polymethylmethacrylate  was  rotated  in  the  illuminating  beam  at  on£* 
revolution  per  minute.  Figure  6  is  a  segment  of  a  dual-channel  chart  recorder  trace  of 
the  output  of  the  two  photomultipliers  recording  these  reflexions  The  upper  trace  i^ 
PMTl  and  the  lower  is  PMT2.  Examination  of  this  figure  shows  that  the  amplitude 
and  shapes  of  pulses  from  different  particles  can  vary  substantially,  and  frequently  the 
individual  pulses  comprising  a  pair  differ  from  each  other. 

As  3hown  in  appendix  A,  the  result  of  these  variations  is  an  uncertainty  in  thf\; 
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measurement  time  between  the  pulses  detected  using  a  threshold  detection  criterion. 
This  ‘  pulse  variation  ’  error  can  be  minimized  by  maximizing  the  separation  in  time  of 
the  two  pulses  relative  to  their  width,  that  is  by  maximizing  the  rotation  angle  a  and 
minimizing  the  diffraction  width  /?,  or  maximizing  the  ratio  a/ fi.  A  narrow  pulse  width 
is  achieved  by  the  use  of  a  short  wavelength  of  light,  for  which  a  He-Cd  laser  (Liconix 
Model  401,  A  —  4431  A)  was  chosen.  The  optimum  rotation  angle  is  difficult  to 
determine  because  all  of  the  deleterious  variations  become  more  pronounced  as  the 
rotation  angle  increases  and  disappear  at  zero  rotation  angle.  Thus  there  is  a  trade-off 
between  maximum  time  separation  and  minimum  shape  variation.  In  practice  choosing 
a//?  ~  1  proved  to  be  satisfactory.  The  measured  elapsed  time  distribution  for  a 
constant  rotation  rate  of  the  dummy  specimen  is  plotted  in  figure  7.  This  distribution 
eventually  determines  the  vorticity  measurement  resolutic  i  to  be  discussed  in  §  3. 

(6)  Electronics 

A  block  diagram  of  the  data  acquisition  system  is  shown  in  figure  8.  Threshold 
detectors  are  employed  to  determine  the  time  of  occurrence  of  the  pulses.  A  pulse  from 
PMT 1 ,  after  amplification,  triggers  a  threshold  detector,  starting  a  clock  which  is  accur¬ 
ate  to  1  pa.  A  similar  event  in  the  second  channel  stops  the  clock  and  disables  the  thres¬ 
hold  detectors.  The  elapsed  time  is  read  by  an  on-line  computer  which  then  resets  the 
clock  and  threshold  detectors,  and  subsequently  operates  on  the  data.  To  suppress  low- 
vorticity  noise,  discussed  in  §  3  below,  the  system  automatically  resets  after  a  pre¬ 
determined  elapsed  time  if  no  pulse  is  presented  to  the  off  detector.  The  data  are  con¬ 
verted  to  vorticity  in  the  computer  and  displayed  on  a  graphics  terminal.  Note  that 
measurements  of  vorticity  of  only  one  algebraic  sign  are  recorded  with  this  system  as 
described.  Both  signs  could,  however,  be  simultaneously  recorded  with  duplicate 
electronics  in  which  the  on  and  off  channels  are  interchanged. 


3.  Analysis  of  system  characteristics 

The  precision  and  accuracy  of  vorticity  measurement,  data  acquisition  rate, 
temporal  resolution,  spatial  resolution  associated  with  the  fluid  volume  sampled,  and 
the  background  noise  spectrum  characterize  the  vorticity  measurement  system.  These 
properties  are  closely  interrelated  and  may  be  optimized  selectively  to  suit  various 
applications. 

A  special  requirement  of  the  optimization  procedure  is  the  selection  of  conditions  that 
facilitate  convenient  separation  of  the  measured  vorticity  spectrum  from  the  charac¬ 
teristic  background  noise  of  the  VOP  system,  which  is  unique  in  that  the  noise  spectrum 
depends  on  the  vorticity  (signal)  spectrum  as  well  as  instrumental  parameters.  Below, 
in  §3.1,  expressions  describing  the  data  rate  and  spatial  and  temporal  resolutions  are 
derived,  the  background  noise  is  discussed  in  detail,  and  a  simple  demonstration  of  the 
optimization  technique  is  presented. 

Section  3.2  describes  the  experimental  resolution  which  can  be  expected  after 
optimization.  A  system  resolution  function  combining  the  effects  of  flow  inhomo¬ 
geneities  and  optical  variations  is  developed  to  facilitate  analysis  of  experimental 
results. 
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3.1.  Data  rate  and  noise  analysis 

(a)  Data  rate.  Assume  a  sampled  volume  of  cross-sectional  area  A  and  uniform  length  «,yj 
a  mirror  concentration  n,  flow  at  mean  velocity  v,  and  a  vorticity  vector  («,,  0, 0). 
The  rate  at  which  mirrors  enter  and  pass  through  the  sampled  volume  is  nAv,  of  which*, 
a  small  fraction,/,  will  present  the  proper  vertical  orientation  to  reflect  light  in  the% 
direction  of  slits  1  and  2,  which  are  positioned  to  measure  the  x  component  of  vorticity. 
Thus  the  effective  mirror  passage  rate  is  Rp  -  nAvf.  The  particles  rotate  through  an, 
angle  tf>  =  ^(s/v)<ol  while  they  transit  the  sample  volume.  Those  that  rotate  through V- 
the  appropriate  azimuth  defined  by  slit  1  and  continue  through  an  additional  angle"*- 
2a  before  leaving  the  sampled  volume  can  generate  vorticity  data. 

The  temporal  resolution  associated  with  an  individual  datum  is  the  time  for.£ 
rotation  through  2a,  that  is  rx  —  4a/w1.  Ideally  a  -*  0  would  yield  instantaneous  data, r-*' 
In  this  limit  the  probability  of  a  mirror  entering  the  sampled  volume  with  an  initial 
azimuth  suitable  for  data  production  is  <j>/2rt,  yielding  an  ideal  limiting  data  rates'^ 
Rt  =  2Rp(j>/2i:  (see  appendix  B),  or 

Ri  =  fnAs(ii1/2n  —  fNw1j2n,  (3.1) 

where  N  is  the  average  number  of  mirrors  within  the  sampled  volume.  -R,  is  also  the 
limiting  data  rate  generated  by  a  mirror  density  n  as  the  flow  velocity  vanishes,  or 
</>-*■  co.  The  mirrors  rotate  at  a  frequency  u)1/tn  with  a  fraction  /oriented  to  produce 
two  data  per  rotation,  due  to  the  twofold  symmetry  of  the  mirrors. 

To  measure  a  rotation  time  the  detector  slits  1  and  2  must  be  separated ;  hence  a  is  v*- 
necessarily  non-zero.  Consequently,  at  non-zero  flow  velocity,  some  mirrors  enter  the 
sampled  volume  with  an  azimuth  which  will  not  allow  data  generation  during  their  •/ 
residence  time.  A  mirror  may  then  generate  pulses  which  either  trigger  the  on  detector  ‘.'j 
and  the  off  detector  in  sequence,  producing  a  ‘  valid  ’  datum,  trigger  only  one  detector, 
or  trigger  no  detector  before  leaving  the  sampled  volume.  Appendix  B  enumerates  , 
these  possibilities  and  shows  that  mirrors  having  an  initial  azimuth  suitable  for  valid  ’-V 
data  production  are  detected  at  a  rate 

Ra  =  fnAv(<f>  -  2a) /n,  (<j>~2a)>0A  < 

Ra  —  0,  (0  —  2a)  <  0.)  (3,  )  - 


The  rate  at  which  the  individual  detectors  are  triggered,  however,  remains  equal  to 
Rit  independent  of  a  and  v.  The  difference,  Rr  =  Ri-  Ra,  defines  the  rate  at  which  pulses 
appear  in  each  detection  channel  that  have  no  correlated  pulses  (those  due  to  the  same 
reflexion)  in  the  other  channel.  These  ‘random  ’  pulses  result  in  a  further  reduction  of 
the  measured  data  rate  by  occasionally  appearing  between  the  two  correlated  or 
‘valid’  pulses  comprising  the  ‘valid  on-off  pair’  required  to  generate  a  valid  datum. 
Random  pulses  occurring  sequentially  in  the  two  channels  may  also  produce 1  false  ’ 
data  values,  resulting  in  ‘  background  noise 

The  measured  valid  data  rate,  Rm,  is  determined  by  enumeration  of  the  possible 
sequences  of  events  involving  valid  and  random  pulses.  Five  such  sequences  can  occur: 

(1)  A  valid  on-off  pair  is  uninterrupted  and  produces  a  valid  datum. 

(2)  A  valid  pulse  starts  the  clock  and  a  random  pulse  stops  it,  generating  noise. 

(3)  A  random  pulse  starts  the  clock  and  a  valid  off  pulse  stops  it,  generating  noise 
and  disqualifying  the  valid  pair  from  data  production. 

(4)  Random  pulses  both  start  and  stop  the  clock,  generating  noise. 
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(5)  Two  (or  more)  valid  on-off  pairs  overlap,  disqualifying  the  second  pairfromdata 
production. 

Sequences  2,  3  and  5  limit  the  measured  data  rate  to  some  value  less  than  Ra.  In 
addition,  sequences  2,  3  and  4  generate  invalid  elapsed-time  measurements,  thus 
creating  a  spectrum  of  false  vorticity  noise. 

The  measured  data  rate  then  is 


Rm  =  Ra-Ri-R3-Ri,  (3.3) 

where  Rj  are  the  rates  at  which  the  enumerated  sequences  j  =  2,  3,  5  occur.  These  rates 
are  calculated  in  appendix  C  with  the  results 


(3.4a) 

(3.46) 

(3.4c) 


(3.5) 


Rt  =  I2a[l  -  exp  ( -  4afl(./<y1)3, 

•ff3  =  RaRf/(Rr  +  Rt), 

Rs  =  -RJ1  “  exp  ( -  4aRJ<o1)]. 

On  inserting  these  results  in  (3.3),  the  measured  data  rate  becomes 

Rm  =  -RJexp  ( ~  +  exp  ( -  4 aRJ  wj  -  (2  Rr  +  Ri)/{Rr  +  )] 

(b)  Noise.  Sequences  2,  3  and  4  above  yield  continuous  distributions  of  values  of  the 
elapsed  time.  Since  these  are  all  Poisson  processes  (appendix  C),  the  probability  density 
of  elapsed  times  is  directly  calculable.  Weighting  the  probability  densities  with  the 
corresponding  noise  generation  rates  yields  the  total  elapsed-time  noise  distribution : 

(RaRrexp(-Rrtn),  tn  <  4a/w1, 
tn  >  4 a/to 


Dn(tn)  =  [RrR\/(Rr  +  *,)]  exp  ( -  Rttn)  +  J*a 


(3.6) 


The  first  term  is  generated  by  sequences  3  and  4,  and  the  second  term  is  by  sequence  2. 
This  distribution  determines  a  fictitious  vorticity  distribution,  or  vorticity  noise 
spectrum,  conveniently  written  as 

B»(un)  =  [4a(i21  —  Ra)/ w*  ]  [exp  ( —  4ai?(/  ton)] 


X 


[R\/(2Ri-Ra) 


,  I Ra  exp  (4aRa/(on)] , 

\0], 


(On  >  <Ov 
<  Vi; 


(3.7) 


This  distribution  has  a  sharp  peak  at  oin(peak)  =  2aRi,  and  a  long  tail  towards  high 
values  of  vorticity.  The  total  noise  data  rate  is 

Rn  =  (R*i  +  RiRa-R*)/(2Ri-Ra)-Raexp[-4oc(Ri-Ra)/<o1].  (3.8) 

(c)  Discussion  and  optimization.  The  measured  data  rate  given  by  (3.5)  is  implicitly 
a  function  of  all  experimental  parameters,  Rm  =  Rm(f,  N,  a,  s,  v,  (ox).  The  VOP  system 
is  optimized  by  maximizing  the  data  rate  for  a  particular  experimental  vorticity  value 
and  spatial  resolution.  The  spatial  resolution  is  limited  by  the  requirement  that  the 
mirrors  have  sufficient  residence  time  in  the  sampled  volume  to  rotate  through  the 
angle  2a  before  translating  out  of  the  sampled  volume,  that  is,  Ra  >  0,  or  by  (3.2) 

Wj/v  >  4 a/s.  (3.9) 

The  data  rate  is  limited  by  overlapping  pulse  pairs  if  the  flow  velocity  is  small,  and 
by  the  requirement  that  the  true  vorticity  peak  be  easily  distinguished  from  the  noise 
spectrum  if  operating  near  the  vorticitv/velocity  threshold  limit  given  by  (3.9).  Since 
both  the  data  rate  and  the  noise  spectrum  depend  on  the  number  of  mirrors  within  the 
sampled  volume,  optimization  generally  imposes  a  condition  on  N.  Examination  of 
equations  (3.5)  and  (3.7)  in  the  limits  v  ->  Oand  <oJv  4 a/s  illustrate  these  limitations. 
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(i)  Zero  flow  velocity  limit.  Since  Rr  -*•  0  and  Ra  -*•  as  v  -*  0,  (3.5)  shows  that  fj. 


lim  Rm  =  R{  exp  ( —  4ai?,/w1). 

tr-*0 


(3.10) 


This  data  measurement  rate  peaks  at  Rt  =  «u1/4a,  resulting  in  a  maximum  possible^ 
data  rate  of 

RmmAx  ='«i/4ea  =  wi/10-87a.  (3.11)$? 

Optimization  of  the  number  of  mirrors  in  the  sampled  volume  for  this  data  rate'*-" 
requires,  from  (3.1), 

iVopt  =  n/2fea  =  0-58//a.  (3. 12)$ 

r> 

Either  more  or  fewer  mirrors  in  the  sampled  volume  will  reduce  the  data  rate. 

(ii)  Vorticity  threshold  limit.  In  the  vorticity  measurement  threshold  limit.  _ 
wjv  -*■  4 a/a,  the  valid  pulse  pair  rate  Rn~*  0  and  Rn  ->  R{,  so  that  all  of  the  pulse^J 
eventually  become  random,  and  the  vorticity  peak  vanishes  into  the  noise  background.**^ 
Appropriate  approximations  of  (3.5)  and  (3.7)  provide  an  analysis  of  the  approach  to 
this  limit.  V 

The  background  noise  spectrum  at  un  =  w1  here  becomes  < 

Bn(oJx)  ~  2 xRyoj*  =  WfX/nf,  (3.13U 

where  the  assumption  2 fNa/n  ^  1  has  been  invoked  to  approximate  the  exponential 
in  (3,7)  by  unity.  Equation  (3.5)  becomes,  approximately, 


Rm  sc  {fNf  2n)  (Wj/2  -2av/s). 


(3-14  , 


Assuming  that  the  valid  data  peak  at  w1  is  actually  spread  over  a  width  Sw,  the  spectral* 
density  of  the  valid  data  is  Rm/8w  so  that  a  ‘  signal  to  noise  ’  ratio  can  be  defined  as  *** 


rn  -  (/*m/<to>)/-B»K)  =  n(wj2-2 av/s)/afNSo>. 


(3.15$ 


Selection  of  a  signal-to-noise  ratio  rn  as  large  as  desired  specifies  a  maximum  valu^, 
ot  n, 

N  ^  n((oj2-2ocv/s)/afrnS(ij.  (3  16f 

With  the  maximum  value  of  N,  the  noise  peak  occurs  at 

V 

&>n(peak)  =  2aR{  =  (ol(w1/2-2av/s)/rnS(t).  (3. 17) 

If  wx/S(i)  <  1  then  w„(peak)  (ox  and  the  false  peak  is  clearly  separated  from  the  re£$ 
peak,  as  desired.  These  results  are  graphically  illustrated  in  figure  9  where  the  valid  pea'v 
is  represented  by  a  gaussian  ot  width  $u>,  rn  =  20,  N  is  determined  by  (3.15)  and  the 
other  experimental  parameters  are  chosen  to  correspond  to  actual  experiments.  ./ 

(iii)  General  case.  The  entire  discussion  of  data  rate  and  noise  has,  to  this  poin>". 
assumed  a  single  value  (or  narrow  range)  of  the  vorticity  vector,  enabling  this  com¬ 
prehensive  illustration  of  optimization.  However,  real  flow  fields  may  have  a  distribiv 
tion  of  vorticity  vectors  spread  out  in  magnitude  and  containing  components  othc^-, 
than  the  one  which  is  sampled.  In  such  situations  the  values  of  Ri  and  Ra  are  still1 
proportional  to  the  mirror  density  re,  but  quantitative  analysis  becomes  excessively 
tedious.  However,  the  preceding  discussion  suggests  that  the  true  vorticity  spectra;  ( 
can  always  be  separated  from  the  noise  spectrum  by  reducing  the  particle  concentration? 
until  the  peak  of  the  noise  spectrum  is  located  at  a  value  of  vorticity  which  is  much 
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Figure  9.  Theoretical  noise  distribution  near  vorticity/velocity  threshold  value. 


less  than  the  threshold  value  of  measurable  vorticity,  wm)n,  for  the  particular  flow 
situation.  This  peak  is  easily  recognized  but  it  can  be  suppressed  by  automatically 
resetting  the  detectors  after  a  time  interval  tr  has  elapsed  after  an  on  pulse  without  an 
off  pulse.  Choosing  tr  so  that  4a/wmin  <^tr  4  l//?{  allows  a  slight  increase  in  the 
measured  data  rate,  and  a  considerable  improvement  in  the  appearance  of  the  vorticity 
distribution  by  suppressing  the  low- vorticity  noise. 


3.2.  Vorticity  resolution 

Vorticity  measurement  uncertainties  limit  the  vorticity  resolution  of  the  VOP  system 
and  broaden  measurements  of  vorticity  distributions.  Limitations  of  optical  geometry 
and  imperfection  of  the  probe  particle  reflexions  introduce  small  errors  in  all  reflexion 
transit  time  measurements  and  thus  in  vorticity  measurements,  whatever  the  flow 
system.  In  addition  finite  spatial  resolution  in  the  presence  of  a  macroscopic  spatial 
vorticity  gradient  broadens  the  vorticity  distribution.  In  order  to  describe  the 
properties  of  the  VOP  system,  a  vorticity  resolution  function  is  defined.  It  consists  of 
a  convolution  of  an  optical  resolution  function  g0  and  a  function  gr  due  to  vorticity 
gradients.  The  contribution  of  noise  is  ignored  in  the  theoretical  analysis  presented 
here. 


(a)  Optical  resolution 

An  experimental  estimate  of  the  optical  resolution  in  the  absence  of  a  vorticity  grad¬ 
ient  is  obtainable  from  the  simulation  vortex  experiment  described  in  §  2.4  (a)  (ii).  The 
measured  time  delay  distribution  shown  in  figure  7  can  be  inverted  to  obtain  the 
measured  optical  resolution  shown  in  figure  10.  This  experiment  determines  the 
optical  resolution  function  <70{to/«j)  of  the  present  VOP  system  for  a  vorticity  spike 
at  cjv  It  displays  a  smooth  optical  resolution  function  peak  with  a  root-mean-square 
vorticity  resolution  of  about  9  %.  Since  this  experiment  suffers  from  a  small  vorticity/ 
velocity  ratio  away  from  the  axis  of  rotation,  it  is  not  surprising  that  a  broad  low-level 
tail  at  higher  vorticity  appears.  It  is  a  combination  of  the  pulse  variation  error 
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Limits  (valid  for  y  >  0,  d  >  c) 

(d-c)/ 2V'2  >  (y-ye)  >  -(d-c)/2J2 
(d  +  c)/2J2  >  (y-yc)  >  (d-c)/ 2V'2 
-(d  +  c)/2v2  <  (y  —  yt)  <  (d-c)/ 2J2 


V(2)  c-2[y-yc-{d-c)/2^2] 
V(2)  c  +  2[y  —  ye+  (d  —  c)/2y'2] 


Table  1 
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Figure  12.  Theoretical  and  experimental  pulse  detection  rates  as  a  function  of  intensity, 
used  to  determine  the  reflexion  intensity  distribution. 

discussed  in  appendix  A  and  the  false  data  background  discussed  in  the  previous 
section.  Note  that  the  pulse  variation  error,  as  well  as  the  noise,  increases  in  magnitude 
as  <f>  -►  2a  since  the  particles  translate  through  a  larger  segment  of  the  gaussian 
incident  intensity  profile.  This  means  that  the  accuracy  of  the  data  improves  as  <jj/v 
increases,  and  that  the  relative  magnitude  of  the  tail  will  depend  on  oj/v. 

(b)  Vorticity  gradient  resolution 

The  measured  vorticity  distribution  in  a  system  with  non-uniform  vorticity  is  deter¬ 
mined  primarily  by  the  size  and  shape  of  the  sampled  volume,  and  the  vorticity  and 
velocity  gradients.  For  the  optical  geometry  shown  in  figure  1 1  the  sampled  volume  is 
a  rectangular  box  oriented  at  45'’  relative  to  a  flow  channel.  The  length  s  varies  with 
position.  Calling  y  the  co-ordinate  relative  to  the  centre  of  the  channel  and  yc  the 
position  of  the  centre  of  the  sampled  volume  relative  to  the  centre  of  the  How  channel, 
then  the  length  at  position  y,  s(y,  yc,  c,  d)  has  the  values  summarized  in  table  1 ,  where 
c  is  the  width  of  the  incident  beam  and  d  the  width  of  slit  3.  The  local  contribution  to 
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the  data  rate,  now  a  function  of  y,  is  approximated  by  (3.2)  with  the  assumption  that 
the  effects  of  random  pulses  and  overlapping  valid  pulse  pairs  are  negligible.  Thus 

Ra  =  RaWuiV’  Vc)>  v(y>  yc)>  s(y>  Vo  c> rf)}]  =  fnAv{<j>  -2a), /tt,  <j>  >  2a,  (3.18) 

where  the  function  is  the  angle  through  which  a  particle  rotates  during  its  transit 
of  the  sampled  volume,  defined  in  §3.1  (a).  A  measurement  of  vorticity  with  the 
sampled  volume  centred  at  point  yc  will  now  result  in  a  vorticity  distribution 
gv(w,  v,c,d),  where  gBx.  Ra.  This  piece  of  the  resolution  function  is  called  the  vorticity  - 
sampling  function. 

(c)  Optically  corrected  sampling  function 

The  non-uniformity  of  the  illuminating  laser  beam  and  the  distribution  of  probe 
particle  reflecting  power  couple  in  a  curious  fashion.  The  illumination  is  non-uniform 
because  the  laser  beam  has  a  gaussian  intensity  profile.  Since  larger  probe  mirrors 
reflect  more  strongly  than  small  mirrors  they  trigger  the  threshold  reflexion  detectors 
at  lower  illumination  intensities  and  thus  are  observable  farther  out  into  the  fringes  of 
the  illuminating  beam. 

Defining  m(I)  dl  as  the  number  of  additional  reflexions  that  will  become  detectable 
as  the  incident  intensity  increases  from  /  to  I  +  dl,  the  total  number  of  mirrors  with 

reflexions  bright  enough  to  be  detectable  isl  m(l)dl  if  the  maximum  incident 

intensity  is  /„.  Defining  the  intensity-dependent  beam  width,  c  =  c(/),  the  numberof 
single  channel  reflexions  detected  per  unit  time  is 


M  =  B  f  c(I)m{I)dI, 
Jo 


(3.19) 


where  B  is  some  constant.  For  a  gaussian  laser  beam  of  width  c0  at  ///„  =  e-2  the  beam 
width  between  points  of  intensity  I  is  c(/)/c0  =  (Jin  (/„//))*  so  that 


M(I0)  =  c0Bj\lln(I0/I)im(I)dI. 


(3.20) 


Fitting  m(I)  to  a  gaussian  function  for  simplicity  yields 


M(I0)  =  c0£  j*^°(J  In  (/0/ 7))lexp  [  —  (///„  — a  )2/262]rf/.  (3.21) 

Figure  12  shows  several  scaled  curves  of  M(I0)vs.I0  with  various  values  of  a  and  6. 
Actual  data,  obtained  by  counting  the  number  of  pulses  triggering  threshold  detector 
1  in  a  fixed  period  of  time  for  various  attenuations  of  the  incident  beam,  most  closely 
fits  the  curve  with  a  =  I,  b  =  0-5.  The  vorticity  gradient  function  gt.  can  now  be 
adjusted  to  g'v(to,  yc,  c0,  d)  to  account  for  the  variations  in  probe  particle  reflecting  power 
by  performing  the  integral 


ft  0 

g’v(u)tyc,c0,d)  =  I  gv[ut,yL.,c(I),d]m(I)dI. 
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Vorticity  (s'1 ) 

Figure  13.  Superposition  of  vorticity  resolution  function  and  measured  distribution.  Flow 
parameters:  (a)  w  =  0-22  cm,  yc  =  0-157  cm,  c0  =  0-1  cm,  d  —  0-1  cm,  v0  —  8-65  cm  s-1, 
(6)  w  =  0-8  cm,  yc  —  0-74  cm,  ca  =  0-1  cm,  d  =  0-1  cm,  v0  =  10-3  cm  s_1. 

(d)  Vorticity  resolution  function 

The  complete  vorticity  resolution  function  p(w,  yc,  c0,  d)  is  the  convolution  of  the 
corrected  sampling  function  g'v  with  the  measured  optical  resolution  function  g0. 
Thus  the  vorticity  resolution  function  is 

g(w,  yc,  c0,  d)  =  f  dojag'(oja,yc,cofd)g0(oj/(oa) 

0 

dM„g„Uo/(oa)  f  dl  </, //,.,  c(/),  d)exu  [  —  2(///0  -  l)2]. 
u  J  u 
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(e)  Illustrative  calculation  of  vorticity  .  ‘snl  t ion  function 
The  simple  example  of  laminar  two-dimensional,  Poiaeuille  flow  (Batchelor  1967, 
§4.2)  illustrates  the  calculation  of  the  vorticity  resolution  function.  The  velocity 
and  vorticity  profiles  are  for  this  case 

v(y)  =  »0  (l-y*/w1),  u>(y)  =  2v9y/wi,  (3.24),  (3.25) 

where  v9  is  the  velocity  at  the  centre  of  the  channel  and  w  is  the  channel  width.  The 
angle  of  particle  rotation  during  transit  of  the  sampled  volume  is 

f>(y,  yc>  c,d)  =  s(y,  yc,  c,  d) y/(u> 8  -  y*).  (3.26) 

The  local  data  rate  within  the  sampled  volume  is 

Ra  =  (fnv0hdy/n)(l-y2/wi)[sy/{u;*-y2)-2a],  (3.27) 

where  h  dy  is  the  cross-sectional  area.  Treating  y  as  a  function  of  vorticity ,  the  vorticity - 
sampling  function  is 

gv(aj,yc,c,d)  =  3[y(<u),y<.,c,d]y(o))-2a[wi-y2(oj)].  (3.28) 

Now  the  vorticity  resolution  function  g(ui,  yc,  c0,  d)  can  be  calculated  since  all  ol  its 
components  are  determined.  The  bold  lines  in  figures  13  (a,  b)  show  the  results  of 
numerical  calculations  of  g(u>(v0,  w),  yc,  c0,  d)  from  (3.22)  for  the  present  VOP  system 
with  the  parameters  v0,  w,  y^  c0  and  d  indicated  in  the  figure  caption.  For  this  high- 
symmetry  geometry  (c0  ~d)  the  vorticity  distribution  is  nearly  symmetric,  except  for 
the  noise  tail.  Therefore  several  convenient  measures  of  the  position  of  the  vorticity 
peak  nearly  coincide,  namely  the  putative  vorticity  at  the  centre  of  the  measuring 
volume  w(yc),  the  peak  vorticity  u>p,  and  the  mean  of  the  half-maximum  points  ojh. 

4.  VOP  system  test 

4.1.  Test  system 

Laminar,  two-dimensional  Poiseuille  flow  provides  a  suitable,  stringent  test  of  the 
reliability  of  the  new  vorticity  measurement  system.  A  Plexiglas  channel,  5- 1  cm  high, 
1-3  cm  wide,  and  50  cm  long  was  constructed  for  this  purpose.  An  extra  piece  of 
Plexiglas,  5-1  cm  high  and  25  cm  long,  can  be  inserted  to  reduce  the  channel  width  to 
a  nominal  0-4  cm.  A  gravity-fed,  recirculating  pumping  system  drives  the  fluid,  dibutyl 
phthalate,  through  the  free-surfaced  channel.  Screens  placed  in  the  entrance  region 
smooth  out  gross  flow  irregularities.  Flow  velocity  maxima  were  determined  by 
measuring  the  transit  times  of  visible  particles  between  fiducial  marks.  The  flow  system 
has  provisions  to  vary  all  of  the  parameters  entering  into  the  resolution  function  of 
§  3,  to  facilitate  its  use  for  testing  and  calibrating  the  VOP  system  and  determining  its 
vorticity  resolutions.  A  cylindrical  lens  placed  in  the  path  of  the  incident  beam  can 
change  the  otherwise  cylindrical  sampled  volume  into  a  rectangular  parallelepiped 
Flow  visualization  was  used  to  determine  the  areas  within  the  channel  where  the 
flow  is  sufficiently  well  developed  for  good  measurements.  Suspended  guanine  platelets 
(Kalliroscope  Corporation  product  AQ- 1000)  are  aligned  by  the  flow  in  the  direction  of 
the  maximum  shear  rate,  making  boundary  layers  and  flow  inhomogeneities  clearly 
visible.  The  velocity  and  vorticity  profiles  become  fully  developed  slightly  downstream 
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Figure  14.  Measured  vorticity  profiles  in  Foiseuille  flow.  Calculated  speeds:  (a)  2-9  cm  s~‘ ; 
(b)  4-8  cm  s_1 ;  (c)  8-6  cm  s-1;  (d)  13-0  cm  s-1;  (e)  10-3  cm  S' ;. 

of  the  point  at  which  the  boundary  layers  from  either  side  of  the  channel  meet.  In  the 
0-4  cm  wide  channel,  observations  show  that  this  occurs  at  positions  at  least  20  cm 
downstream  of  the  entrance.  In  the  1-3  cm  wide  channel  it  was  not  clear  that  the 
Poiseuille  flow  developed  fully  for  velocities  above  3  cm  s~l.  Effects  of  the  bottom  of 
the  channel  appeared  to  be  negligible  above  2-5  cm.  Velocities  at  the  centre  of  the 
channel  range  from  about  1  to  15  cm  s-1. 

4.2.  Vorticity  distributions 

Figure  13(a)  compares  an  anticipated  vorticity  distribution  calculated  in  §  3  with  the 
corresponding  measured  distribution  using  the  same  flow  parameters  in  the  narrow 
channel.  The  agreement  is  quite  good,  especially  near  the  half-maximum  points.  The 
noise  and  pulse  variation  errors  for  this  flow  are  apparently  less  significant  than  in  the 
simulation  vortex  experiment  of  §  2.4(a)  (ii),  as  indicated  by  the  differences  in  the  tail 
heights.  Figure  13  (6)  shows  a  similar  comparison  of  some  wide-channel  data  taken  near 
the  channel  wall.  Note  that  the  wide  channel  has  a  shallower  vorticity  gradient  than 
the  narrow  channei,  with  a  correspondingly  narrower  sampling  function  as  expected. 
It  also  has  a  smaller  vorticity/velocity  ratio,  which  is  responsible  for  the  increased 
experimental  tail  amplitude  and  the  small  noise  peak  at  low  vorticity  values.  This 
distribution  can  be  compared  to  the  theoretical  noise  distribution  of  figure  9.  The 
automatic  reset  circuit  is  responsible  for  the  sharp  cut-off  at  the  low  end  of  the  experi¬ 
mental  distribution.  The  experimental  wide-channel  data  shows  a  curious  secondary 
peak  centred  at  approximately  twice  the  value  of  the  primary  peak  which  appears  only 
in  observations  close  to  the  wall.  It  may  be  the  result  of  some  sort  of  flow  fluctuations. 
The  Reynolds  number  ( vawfv )  for  this  flow  is  about  100. 

A  numerical  device  which  is  useful  to  compare  different  vorticity  distributions  is 
demonstrated  by  the  rectangle  drawn  in  figure  13(a).  The  width  and  position  of  this 
rectangular  construction  are  determined  to  yield  the  second  and  third  moments  of  the 
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experimental  distribution  around  its  average.  The  bottom  10  %  of  the  distribution  is 
ignored  in  calculating  the  moments  to  eliminate  the  contributions  of  the  various  errors, 

(a)  Vorti  ity  profiles 

Several  sets  of  vorticity  measurements  are  plotted  as  functions  of  measurement 
position  yc  in  figure  14.  They  show,  within  error  limits,  the  anticipated  linear  vorticity 
profiles  for  these  Poiseuille  flows,  (3.25).  Each  vorticity  point  is  plotted  as  a  rectangle 
with  a  width  representative  of  the  positional  uncertainty  and  a  height  equal  to  the 
difference  between  two  measures  of  vorticity,  namely  the  peak  vorticity  wp  and  the 
average  of  the  half-maxima  wA.  In  accord  with  §  3.2  (d),  wp  and  wA  should  very  nearly 
coincide  with  the  vorticity  at  the  centre  of  the  sampled  volume,  w(yc).  The  experi¬ 
mental  values  of  wA  are  repeatable  to  within  a  few  per  cent  and  are  thus  good  measures 
of  oj(yc).  The  r.m.s.  value  of  the  differences  wp  —  wA  provide  a  practical  measure  of  the 
uncertainties  of  vorticity  measurements.  For  these  data  the  relative  uncertainty  of 
vorticity  is  about  8%. 

The  velocities  indicated  in  figure  14  are  calculated  from  the  slopes  of  the  straight 
lines  determined  by  the  vorticity  measurements,  since  this  measure  is  more  precise 
than  the  flow  velocity  estimates.  The  calculated  velocities  do  agree  with  the  approxi¬ 
mate  measured  velocities. 

In  the  wide  channel  the  flow  has  not  reached  a  fully  developed  Poiseuille  parabolic 
velocity  profile.  Boundary-layer  theory  (Landau  &  Lifshitz  1959,  ch.  4)  shows  that  the 
boundary  layer  20  cm  downstream  from  the  edge  of  a  flat  plate  will  have  a  width  of 
about  0-5  cm  when  the  flow  velocity  is  1 0  cm  s"1  and  kinematic  viscosity  is  0- 1 5  cm2  s~' . 
Thus  there  should  be  a  high- vorticity  boundary  layer  with  small  vorticity  gradient  near 
the  wall  that  drops  rapidly  to  zero  near  the  edge  of  the  boundary  layer.  The  wide- 
channel  vorticity  profile,  line  ( e )  in  fig.  14.  appears  to  have  a  fully  developed  flow 
profile  near  the  wall  for  a  channel  width  of  1  -6  cm  at  the  correct  centre-line  velocity  of 
10- 3  cm  s_1.  However  the  actual  channel  width  is  only  1-3  cm  and  the  vorticity  drops 
sharply  below  the  detectable  minimum  beyond  0-25 cm  from  the  wall.  These  results 
are  consistent  with  the  above  calculations  and  with  the  flow  visualization  observations 
that  a  distinct  boundary  layer  remains. 

(b)  Distribution  widths 

Poiseuille  flow  allows  a  very  simple  comparison  to  be  made  between  the  expected 
vorticity  resolution  and  the  actual  width  of  the  distribution.  According  to  equation 
(3.25),  if  a  width,  Ay,  of  the  vorticity  profile  is  uniformly  sampled  (s  =  constant),  then 
the  corresponding  width  Aw  of  the  vorticity  distribution  is 


a 


A  w  =  (2r0/ic2)  A  y. 


(4.1) 


.v 


In  the  present  experiment  (4.1)  should  provide  a  reasonable  approximation  even 
though  3  varies.  Figure  15  shows  a  plot  of  mean  values  of  Aw  for  each  of  the  lines  in  ‘A> 
figure  14  vs.  2 v0/w2  and  a  straight  line  fit  for  comparison  with  A y  determined  by  (4.1). 
Here  Aw  is  chosen  as  the  width  of  the  rectangular  construction  discussed  above  The 
resulting  value  for  Ay,  ~  0- 1  cm,  agrees  with  the  value  determined  by  the  apparatus  » 
optics. 
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Figure  15.  Experimental  vorticity  distribution  widths  as  a  function 
of  flow  parameter.  Slope  =  0-09  cm. 


5.  Summary 

A  system  for  vorticity  measurement  by  direct  optical  probe  has  been  described  and 
tested  under  the  stringent  conditions  imposed  by  laminar  flows.  The  accuracy,  vorticity 
resolution,  data  rate  and  spatial  resolution  of  the  VOP  and  the  experimental  para¬ 
meters  on  which  these  properties  depend  have  been  analysed  and  demonstrated  by  the 
measurements  with  a  rather  primitive  apparatus. 

There  are  some  fundamental  interdependent  limits  on  the  VOP  system  capabilities. 
The  drift  velocity/ vorticity  ratio  sets  a  lower  limit  on  the  spatial  resolution, 
Smin  =  4 av/uj,  since  the  residence  time  of  a  probe  particle  within  the  sampled  volume 
must  exceed  the  vorticity  detection  time.  An  upper  limit  on  the  number  of  particles 
within  the  sampled  volume  with  zero  flow  velocity,  Nov t  =  0-58 //a,  is  determined  by 
maximizing  the  data  rate  while  requiring  that  reflexions  from  multiple  particles  have 
a  small  probability  of  overlapping.  The  factor  /  represents  the  geometrical  data  col¬ 
lection  efficiency,  which  may  be  improved  by  more  elaborate  optical  design.  The 
maximum  data  acquisition  rate,i?m  msx  =  w/10-87a,  is  achieved  in  a  zero  flow  velocity 
situation  if  the  number  density  is  determined  iV0pt.  As  the  flow  velocity  increases  from 
zero  the  maximum  data  rate  decreases  and  background  noise  increases.  The  signal  to 
noise  ratio,  rn  x  rr(  Jw  —  2 ocv/s)/a.fN8u>  can  be  made  as  large  as  desired  for  any  value  of 
vorticity  greater  than  the  limiting  value  of  to  =  4 ccv/s  by  selecting  a  sufficiently  small 
mirror  number  density,  thus  separating  the  valid  data  from  the  background.  Table  2 
summarizes  this  discussion  of  limiting-system  capabilities  and  gives  the  parameters  of 
the  present  realization  of  the  VOP  system. 

The  discussion  of  the  background  noise  reveals  one  noteworthy  point.  As  x  -*■  0  all 
•cjise  disappears,  the  particle  number  density  is  limited  only  by  the  size  of  the  spherical 


VA-'.vv 


sew 


«0 WWWR 


WWW 


VV  VVWVVA'V'V«V,' 


I 


196 


M.  B.  Frisk  and  W.  W.  Webb 


Parameter 

Typical  value 

HV  i 

(A)  Spatial  resolution 

jj 

=  4av/ti> 

a  =  0-05 

1 200  n  (v / u)  =  0-1  cm) 

■y  ' 

8mi°  \900/i(v/oj  =  0-5  cm) 

9.  } 

g 

(B)  Date  rate 

(1)  Maximum 

=  to/10-Sla 

=  l-84w 

[ 

(2)  Actual 

>>  ' 

(a)  Zero  flow  velocity 

£  : 

B,  =  ioj  (fN/n) 

V  =  10 
/  =  0-04 

V  .  i 

f 

II 

05 

V7 

(6)  Near  io/v  threshold 

Ra  =  (i«j-2at>/«)  (JN/n) 

s 

V 

(C)  Signal-to-noise  ratio 

\ 

r„  =  —  2at,v/s)/afNS<i) 

o)  =  22  s-1 
v  =  10-3  cm  s_1 

So)  =  10  8~ 1 

r„  =  22 

{D)  Vorticity  resolution 

i  j 

V  3 

Present  practical  limit 

~  9% 

Tabus  2.  Present  VOP  system  parameters. 

-  v»  1 

probe  particles,  and  the  data  rate  is  bounded  by  only  this  density  limit.  In  the  present  r^>  \ 

VOP  system  a  must  be  non-zero  to  deal  with  variations  of  the  size  and  shape  of  the  'j 

reflexions.  The  rotation  angle  2a  is  coupled  to  the  divergence  angle  of  the  micro-mirror  j 

reflexions  by  the  requirement  a//?  >  1.  However,  if  particles  with  completely  uniform 
mirrors,  preferably  perfect  disks,  could  be  produced,  and  if  the  illuminating  beam  !"*•  jj 
were  converted  from  a  gaussian  to  a  1  top-hat  ’  intensity  profile,  the  vorticity  could  be  i| 

determined  by  measuring  the  width  of  a  single  pulse  peak  as  it  transits  a  slit.  The  r  \ 

minimum  value  of  a  could  then  be  determined  by  the  temporal  resolution  of  the  ’£•  l 

electronics,  and  would  be  very  close  to  zero.  Since  the  vorticity  resolution,  discussed  *  > 

at  length  in  §  3,  is  also  coupled  to  the  optical  parameters,  the  illuminating-beam  , 

intensity  profile,  and  the  particle  uniformity,  it  would  also  be  substantially  improved  •/  ] 

by  better  particles.  Substantial  improvement  of  all  of  the  operational  parameters  -N  •] 

can  be  expected  in  future  VOP  systems.  Nevertheless  the  need  for  refractive  index  J 

matching  of  fluid  and  particles  seems  inescapable.  jj 

There  appear  to  be  no  substantial  obstacles  to  measurement  of  vorticity  spectra  and 
spatial  correlations  even  in  turbulence  although  various  potential  capabilities  of  the  1 
system,  particularly  those  applicable  to  turbulence,  have  not  yet  been  fully  established  ’,■%*  jj 
All  three  components  of  vorticity  should  be  simultaneously  measurable  by  straight-  Jj 

forward  elaboration  of  the  present  VOP  system .  Some  special  properties  of  the  vorticity  <1 

field  are  readily  accessible  since  positive  and  negative  signs  of  vorticity  components  are  iv  5 

automatically  separated.  The  VOP  can  readily  average  vorticity  distributions  over 
large  flow  volumes.  3 
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We  have  adapted  the  VOP  system  to  vorticity  measurements  in  mixing  layers, 
wakes  and  turbulent  boundary  layers  and  are  proceeding  to  investigate  grid  turbulence 
and  coherent -structure  phenomena  of  the  type  analysed  in  recent  experimental  and 
theoretical  works  (Roshko  1976;  Kovasznay  1978;  Aref  &  Siggia  1980). 
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Appendix  A 

Consider  the  measurement  of  the  time  interval  between  two  roughly  gaussian- 
shaped  pulses  which  may  have  truncated  or  very  noisy  tops.  The  noise  prevents  the 
use  of  a  peak  detector  as  an  accurate  indicator  of  the  occurrence  of  an  event.  A 
threshold  detector  may  provide  more  reliable  data,  but  introduces  errors  if  the 
pulse  amplitudes  and  widths  are  not  constant. 

Assume  the  pulses  have  Gaussian  shape  but  different  width  and  height,  with  their 
peaks  separated  by  a  true  time  rp.  The  pulses  are  described  by  the  equations 

Sx(t)  =  ^4 !  exp  [  —  <2/2cr«],  (A  1) 

S2(t)  =  A2exp[(t-rp)2/2(r2].  (A2) 

If  the  threshold  detectors  trigger  at  an  amplitude  a0  the  measured  elapsed  time  is 
t  =  t2-t1  =  Tp  +  o-1(21n  AjI- <7-55(2  In  .4.;)*.  The  fractional  error  is 

At  t-tp  _  cr1(21nA1)I  — <r2(21n (A  3) 

^*P  ^  p  ^  p 

The  error  is  minimized  by  maximizing  rp  and  minimizing  crl  and  <r2.  If  crx  x  cr2  and 
Ax  =  CA2,  where  C  is  a  constant,  then  At  x  <72>/(2)  [(In  C  +  In  A2)b  —  (In  A2)l],  which 
approaches  zero  as  A2  increases,  suggesting  that  large  pulse  amplitudes  are  also 
desirable. 

Appendix  B 

Here  are  derived  the  rates  of  registry  of  valid  vorticity  data  and  of  false  data  due  to 
velocity-induced  rapid  transits  of  the  sampled  volume,  and  the  resulting  false  vorticity 
distribution. 

Defining  a  co-ordinate  system  for  mirror  orientation  in  which  the  azimuthal  direction 
to  the  point  half-way  between  slits  1  and  2  is  taken  as  \n,  a  mirror  enters  the  sampled 
volume  with  azimuth  d0  and  exits  at  azimuth  Q  =  da  +  <j),  where  0  =  ^(s/v)co.  To  register 
a  valid  vorticity  datum  it  must  pass  through  azimuths  \n-  cc  and  then  \n  +  a  to  trigger 
PMTl  (on)  and  PMT2  (off)  in  sequence  (see  figure  1 1).  While  traversing  the  sampled 
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ft  <  2a 

2a  <  <f>  <  2n  —  2a 
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\ix  —  (j>  —  2a 

2n 

t 

K 

t  When  (ft  >  2n  +  2a  substitute  <ft  —  2jn  for  (ft,  where  j  is  the  integer  indicating  the  number  of 
complete  revolutions  that  have  occurred :  j  type-a  events  and  one  additional  event  must  occur 
for  each  particle. 

t  Note  that  Pa+e+Pt,+'  + P„+4  =  (<f>  +  2a)/2ir. 

Table  3 


volume  a  single  rotating  mirror  with  <j>  <  2 jt,  referred  to  below  as  an  event,  can  trigger 
the  detectors  in  several  possible  sequences,  generating  the  following  five  distinct  cases. 

Event  Result 

(а)  0O  <  \n  —  a,  6  >  \n  +  a  good  data 

(б)  0O  >  \n-a.,  0  >  \n  + a  off  only 

(c)  #0  <  \i r  —  a,  0<\n  +  a  on  only 

( d )  0  <\n— a  no  detection 

(e)  0O  <  fa T  +  a,  0  >  in -a  off  followed  by  on 

If  (j>  >  2n  then  each  mirror  can  generate  more  than  one  event.  Table  3  shows  the 
probability  of  occurrences  of  each  type  of  event  generated  by  an  individual  mirror. 
Px  denotes  the  probability  of  only  a  type  x  event,  while  Px¥v  denotes  the  joint  proba¬ 
bility  of  a  type  x  event  followed  by  a  type  y  event. 

The  average  rate  at  which  type  x  events  occur  is  equal  to  twice  the  mean  rate 
(Rp  =  fNAv)  at  which  mirrors  with  proper  (detectable)  vertical  orientation  flow 
through  the  sample  volume  times  the  event  probability  Px, 

Rx  =  2RpPx,  (Bl) 

except  that  type  a  events  (valid  data  generators)  occur  at  the  rate  Ra  =  2 (lc  +  Pa)Rp, 
where  the  integer  k  is  the  number  of  complete  rotations.  The  factor  of  2  is  a  result 
of  the  rotational  symmetry  of  the  mirrors.  Of  interest  are  the  observations  that,  for  all 
0  >  2a,  Ra  =  Rp(<f>  -  2a) /2n,  and  the  rates  at  which  all  on  (type  a  or  c  or  e)  and  all  off 
(type  a  or  b  or  e )  events  occur  equal,  not  surprisingly,  the  ideal  data  rate,  or 

Ron  =  R0ff  —  Ri  =  fy  io/2n. 
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Appendix  C 

The  noise  generation  rates  due  to  the  various  sequences  of  presentations  of  random 
and  valid  pulses  are  calculated  here.  Recall  that  pulses  presented  to  the  off  detector 
before  the  clock  has  been  started  pass  undetected.  Similarly,  pulses  presented  to  the  on 
detector  while  the  clock  is  running  also  pass  undetected. 

Consider  those  sequences  wherein  random  pulses  presented  to  the  on  detector  with 
rate  Rr  start  the  clock,  and  either  random  or  valid  pulses  presented  to  the  off  detector 
at  rate  R{  stop  it.  The  rate  Rf  of  occurrence  of  false  data  generating  pulse  pairs  in  this 
case  is  the  rate  of  random  pulse  presentations  to  the  on  detector,  Rr,  multiplied 
by  the  probability  of  detection  of  these  pulses,  1/(1  +  Rr/Ri),  or 

Rf  =  Rr  R(/ (R(  +  Rt ) .  (Cl) 

This  case  can  be  divided  into  two  sub-cases  where  the  off  detector  is  triggered  hv  a 
valid  pulse,  or  it  is  triggered  by  a  random  pulse.  The  first  has  a  probability  Ra/Ri  of 
occurrence,  and  the  second  has  probability  Rr/Ri- Hence,  the  rate  of  random-onj  +  valid- 
offs  (sequence  3)  is 

R^{RJRi)Rf  =  RTRj{RT  +  Ri)  (C  2) 

and  the  rate  of  random-ww-random-o/fs  (sequence  4)  is 

Rx  =  m/M  +  Ri).  (C  3) 

Random  events  occurring  at  an  average  rate  R  describe  a  Poisson  process  (Lindgren 
1976).  Thus,  starting  at  an  on  detection,  the  normalized  probability  density  of  elapsed 
times  until  a  random  off  detection  is 

P(R,  t)  =  R  exp  ( -  Rt),  (C  4) 

where  the  off  detections  occur  at  rate  R.  For  the  cases  considered  above  the  off  rate  is, 
from  (3.1),  Rit  so  R  =  Rt.  The  resulting  elapsed  time  distribution,  weighted  by  the 
false  data  rate,  is  D{t)  =  RfP(Rit  t). 

Now  consider  sequence  2  where  a  valid  pulse  pair  starts  the  clock  and  a  random  pulse 
stops  it.  There  is  an  additional  restriction:  the  random  pulses  must  occur  in  the  time 
interval  [0,  ia/cjj].  Thus  the  probability  of  this  sequence  occurring  is 

r  *  a/w, 

P2  =  J  P{Rr>  t)dt  =  1  -exp  (  -  iaRr/ojl).  (C5) 

Since  valid  pulse  pairs  occur  at  rate  Ra,  the  rate  at  which  this  source  of  false  data  is 

generated  i3 

Rt  =  RaP2  =  -Sat1  -  exp  ( -  4aRr/(j1)].  (C  6) 

Finally,  the  probability  of  two  valid  pulse  pairs  overlapping  (sequence  5)  is 

P5  =  J^  P(Ra,t)dt  (Cl) 

and  the  rate  at  which  data  is  lost  to  this  sequence  of  events  is 

=  Pnp5  =  #„[!  -exp(-4jt/?n/w,)].  (C’8) 
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c  0MEGA3  VERS.  1  DEC  87 

C 

C  THIS  PROGRAM  READS  A  FILE  OF  DIGITIZED  STREAM  LOCATIONS  ON  FILM 

C  COMPOTES  THETA, PHI  AND  ASSOCIATED  COEFFICIENTS  FROM  X.Z  PAIRS 

C  COMPUTES  dTHET A  /  d  t , dPH I  /dt  _ ACCUMULATES  SUMS,  AND  SOLVES 

C  3X3  TO  EXTRACT  W;:,  Wy,  Wz. 

C 

IMPLICIT  REALMS  (A-H,0 -Z) 

REALMS  PHP I 20) , THR (20) 

REmL*3  SD  <3)  ,  SA  v  5 , 5 )  ,  W  <  --■ ) 

REAL*4  XX, ZZ 
CHARACTER -*30  FLNAME 
PI =3. 14159 

C  LOAD  DATA 

C 

NR  I TE ( * ,  1 00  > 

FEAD <*, 101 ) FLNAME 
C 

C  ENTER  V  COORDINATE  >. DISTANCE  TO  FILM  PLANE' 

C 

2  WRITE ' 300) 

READ <*, 301 ) V 
WRITE  (*,  302.'  1 
WRITE ' *, 398) 

SY=1 . DO 

IF ( Y . LT . 0 . DO ) SY=- 1 . DO 

OPEN ( 1 , F I L£=FLNAME , STATUS=  OLD ' > 

DO  3  1=1,20 

READ ■ 1 , 200, END=4,ERR=999> XX, ZZ 
=  X  X 


RF=DSQRT <X*X+Z*Z> 

TT I =RF  DABS iY) 

TTR=DTAN(DASIN( 1 . 5D0*DS IN < DATAN ( TT I  /  t  >  ) 
RN=RF— O. 4 DO* I TTR-TTI) 

X  =  = *RN/RF 
2  =  Z*RN/ RF 

RHP  < I > =S i  *R I / 2 . DO-DATAN ( x / Y  > 

THR ( I > =R 1/2. DO-DATAN < Z.  DSQRT ( X*x  +  Y*Y>  / 
WRITE!*, 399 ) THR ( I ) , PHP ( I ) 

N=  I 

3  CONTINUE 

4  C_OSE  •  UNI  T  =  2:> 

c 

L  INITIALIZE  ARRAYS 

C 

32=0.  DO 
DO  6  1=1,3 
SD  ( 1:  =0  .  DO 
w  U  .i  =  0 .  DO 
DO  5  10  =1,3 
SA  '  1  , 1-  1  =0 .  DO 

5  CONTINUE 

6  CONTINUE 


C 


ENTER  T IMF  INTERVAL  DT 


$ 

$ 


Wr  I  IP'*,  4',"  ' ) 

*- EaD  ■  *  ,  40 1  •  DT 
WR  ;  TE  (  -  .  4, DT 


ft 


U  u  u  u  o  o  f'  u  u  u  u  n  <*  m  u  u 


COMFUTE  DERIVATIVES,  COEFFICIENTS,  AND  ACCUMULATE  SUMS 
DO  7  L=2.N-1 

DTHFDT  = <  THP (L+l ) -THP (L-l ) ) / (2, DO*DT> 

DFHPDT= (PHP (L+l ) -PHP (L-l ) )  / (2. DO*OT> 

CALL  TRANS ( THP ( L ) , PHP ( L ) , SV , A, B, C , D , E , F ) 

32=S2+DTHPDT +DTHPDT +DPHPDT  *DPHPDT 

SD ( 1 ) =SD ( 1 ) +DTHPDT  *A+DPHF'DT *D 
SD (2) =SD (2) +DTHPDT *B+DPHPDT *E 
SD (3) =SD ( 3 ) +DTHPDT*C+DPHPDT*F 

SA ( 1 , 1 ) =SA (1,1) +A*A+D*D 
SA (1 , 2 ) =SA (1,2) +A*B+D*E 
SA ( 1 , 3) =SA (1,3) +A*C+D*F 
SA (2,2) =SA (2, 2) +B*B+E*E 
SA (2, 3) =3A (2, 3) +B*C+E*F 
SA (3,3) =SA (3, 3) +C*C+F*F 
CONTINUE 
S3 1 1 =SA ( 1 , 1) 

SS22=SA (2, 2) 

SS33=SA ( 3,3) 

S3 1 3=SA (1,3) 

SS23=SA (2,3) 

SS12=SA( 1,2) 


...  A  SYMMETRIC  MATRIX 

SA  (3,  2)  =SA  (2,  3) 

SA (3,1> =SA (1,3) 

SA (2, 1 ) =SA (1,2) 

DO  20  M=  1 , 3 

WRITE (*,  44-9)  SD  <M> , SA (M,  1  )  , SA (M, 2)  , SA (M, 3) 

FORMAT (  0  ' , 1 PD 12 . 4 ,  ' , 3 ( 1PD12. 4) ) 

CALL  MATINV (SA. 3. DET) 

DO  30  M= 1 , 3 

WRITE (*, 450) SA (M,  1 )  ,  SA (M, 2) , SA(M, 3> 

COMFUTE  W::,Wy,Wz 

DO  9  1=1,3 
DO  8  J  =  1 , 3 

W  <  I  1  =W  <  I  ;  +SA ( I , J ) *SD <  J ) 

CONTINUE 

continue 

COMPUTE  CHI-SQUARE  AND  UNCETAINTIES 
NU=2*N-~ 

CH  I  D=W  <  1  >  *W  <  1  )  *35  1 1+W 1 2) *W ( 2  >  *SS22+W  ( 3)  *W  >  3  +SS33 

CHI  C  =  2 .  DO*  <  W  (  1  ■  *W  < 3 >  *SS  1 3+W  •  1  )  +W  ( 2 *SS  1  2+W  >  2  >  *W  •.  3  >  +  3S23 

CH I y=-2. DO* ( W ( 1 ) *SD ( 1 > +W i 2 ) *SD (2 ) +W i 3 ) *SD ( 3 ' > 

CHI  2=  « CHID+CHIC+CHI r+S2)  NU 
WR  I  TE  * ,  499  ■  NU.CHI2 

DW 1 =OSQR T i CH 12/ S3  1  1  / 


I 


c 


s 

99 

999 

73 

C 

C 

% 

j 

100 

101 

C\ 

200 

v> 

300 

kb 

30 1 

302 

Xi» 

398 

fcj 

399 

400 

401 

$ 

402 

450 

499 

500 

50 1 

$ 

6  00 

700 

ft 

C 

H 

c 

c 

v  - 

c 

s 

•*» 

c 

:? 

c 

* 

c 

c 

v; 

■  m 

I 

! 

I 


c 

c 

c 

*r 

c 

r 


« 


DW2=DSQR:T «CHI2. 3S22) 
DW3=DSQRT (CHI2/SS33) 

WRITER,  500)  W(l)  ,  W  (2)  ,  W  ( 3 ) 
WRITE >*,501)  DW1.DW2.DW3 
STOP 

WR  I TE  )  *  ,  600  > 

WR I  TE  *  * ,  7 00  ) 

STOP 


FORMAT <  0  ENTER  INPUT  FILE  NAME:  ’,*) 

FORMAT  ( A 1 5 ) 

FORMAT (2F14. 7) 

FORMAT (  0  ENTER  PERPENDICULAR  DISTANCE  TO  FILM  PLANE:  ,4) 

FORMAT (D12. 4) 

FORMAT (1 PD  12. 4) 

FORMAT (  O  THETA'  PHI  ' > 


FORMAT  ^2  < 1PD12. 4)  ) 

FORMAT (  O  ENTER  CHOPPER  PERIOD:  ,4) 
FORMAT (D12. 4) 

FORMAT (1PD12. 4> 


FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
FORMAT ( 
END 


'0  MATRIX  INVERSE  , 
'  0  CH I -SQUARE  PER 

0  Wz  =  . 1PD12. 4, ' 

'  0  +  /—  ,1F‘D12.4, 

OPEN  ERROR  ) 

READ  ERROR ' ) 


< 1PD12. 4) ) 

,12,  DEGREES  OF  FREEDQM=  . 1PD12. 
W:<=  ,  1PD12.  4,  '  Wy=  ‘  ,  1  PD  12.  4) 

' . 1PD12. 4, '  • . 1PD12. 4) 


4) 


TRANSFORM  TO  UNPRIMED  ANGLES,  COMPUTE  COEFFTS. 


SUBROUTINE  TRANS < THP , PHP , SY , A, B, C , D, E , F > 

IMPLICIT  REAL*8  (A-H.O-Z) 

FAC= 1 . DO+DSIN (THP) *DCQS (PHP) 

TH=DACOS (DCOS (THP) / DSQRT < 2 . DO*FAC ) > 

PH=5Y*DACQS  (FAC/DSQRT  (2.  DO*  (FAC-DCOS  (THP)  *DCOS  (THF  ) /2.D0)  ;  > 


STH=DS I N  <  TH ) 

CTH=DCOS ( TH ) 

SF'H=DS  I N  ( PH  ) 

CPH=DCOS (PH) 

STHP-DSQRT < 1 . DO-4. DO*STH*STH*CTH*CTH*CPH*CPH ) 


COEFFICIENTS 

A  =  2  ■  Do*STH*CTH*SFH/STHF' 

B=-2  ■  DO*SPH*CFH*STH*STH/  STHF' 

C=2. DO* i STH*STH*CPH*CPH-CTH*CTH) /STHF 

D=2. DO*STH*STH* ! 1 . DO-2. DO*CTH*CTH*CPH*CPH) / < STHP+STHP > 

E=2. DO*STH*CTH*CPH* < 1 . DO-2. DO*CFH*CFH*STH*STH )  (STHP*STHP> 
F=-2. DO*STH*CTH*SPH* ( 1 . DO+2. DO*CPH*CPH*STH*STH)  STHP*3THF) 

return 

END 
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